Weak decay of hyper nuclei 



W. M. Alberico and G. Garbarino 

Dipartimento di Fisica Teorica, Universita di Torino and INFN, Sezione di Torino, 
1-10125 Torino, Italy 



1. Introduction 

The focus of these lectures is on the weak decay modes of hypernuclei, with special 
attention to A~hypernuclei. The subject involves many fields of modern theoretical and 
experimental physics, from nuclear structure to the fundamental constituents of matter 
and their interactions. The peculiar behaviour of matter containing strange quarks has 
raised in recent years many interesting problems, one of which being the physics of 
hypernuclei. 

Hypernuclear physics was born in 1952, when the first hypernucleus was observed 
through its decays [1]. Since then, it has been characterized by more and more new 
challenging questions and answers. The interest was further raised by the great advances 
made in the last 15-20 years. Moreover, the existence of hypernuclei gives a new dimen- 
sion to the traditional world of nuclei (states with new symmetries, selection rules, etc). 
They represent the first kind of flavoured nuclei, in the direction of other exotic systems 
(charmed nuclei and so on). 

Hyperons (A, E, S, fi) have lifetimes of the order of 10"""^° sec (apart from the 
which decays into A7). They decay weakly, with a mean free path A « ct = 0(10 cm). 
A hypernucleus is a bound system of neutrons, protons and one or more hyperons. We 
will denote with Z a hypernucleus with Z protons, Z neutrons and a hyperon Y . 
In order to describe the structure of these strange nuclei it is crucial the knowledge of the 
elementary hyperon-nucleon {YN) and hyperon-hyperon {YY) interactions. Hyperon 
masses differ remarkably from the nucleonic mass, hence the flavour SU (3) symmetry is 
broken. The amount of this breaking is a fundamental question in order to understand 
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the baryon-baryon interaction in the strange sector. 

Nowadays, the knowledge of hypernuclear phenomena is rather good, but some open 
problems still remain. The study of this field may help in understanding some important 
questions, related to: 

1. some aspects of the baryon-baryon weak interactions; 

2. the YN and YY strong interactions in the = 1/2+ baryon octet; 

3. the possible existence of di-baryon particles; 

4. the renormalization of hyperon and meson properties in the nuclear medium; 

5. the nuclear structure: for instance, the long standing question of the origin of the 
spin orbit interaction and other aspects of the many-body nuclear dynamics; 

6. the role played by quark degrees of freedom, flavour symmetry and chiral models 
in nuclear and hypernuclear phenomena. 

Many of these aspects can be discussed and understood by investigating the hypernuclear 
weak decays. Important related arguments, which here will not be considered or only 
briefly mentioned, can be found in the lectures by A. Gal [2] and T. Nagae [3], while the 
experimental viewpoint on the same subject is presented by H. Outa [4]. 

In these lectures the various weak decay modes of A-hypernuclei are described: in- 
deed in a nucleus the A can decay by emitting a nucleon and a pion {mesonic mode) as it 
happens in free space, but its (weak) interaction with the nucleons opens new channels, 
customarily indicated as non-mesonic decay modes. These are the dominant decay chan- 
nels of medium-heavy nuclei, where, on the contrary, the mesonic decay is disfavoured by 
Pauli blocking effect on the outgoing nucleon. In particular, one can distinguish between 
one body and two-body induced decays, according whether the hyperon interacts with 
a single nucleon or with a pair of correlated nucleons. 

An interesting rule for the amount of isospin violation (A7 = 1/2) is strongly sug- 
gested by the mesonic decay of free A's, whose branching ratios are almost in the pro- 
portion 2 to 1, according whether a -iT~p or 7r°n are emitted. This totally empirical rule 
has been generally adopted in most of the models proposed for the evaluation of the A- 
hypernuclei decay widths: some of the expected consequences, however, seem to require 
additional work and investigation. Indeed, the total non-mesonic (Fnm = ^n+^p (+1^2)) 
and mesonic (Fm = T^o -|- F^- ) decay rates are well explained by several calculations; 
however, for many years the main open problem in the decay of A-hypernuclei has been 
the discrepancy between theoretical and experimental values of the ratio F„/Fp. This 
topic will be discussed at length here, together with the most recent indications toward 
a solution of the puzzle. 

Another interesting and open question concerns the asymmetric non-mesonic decay 
of polarized hypernuclei: strong inconsistencies appear already among data. Also in 
this for the F„/Fp puzzle, one can expect important progress from the present 
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and future improved experiments, which will provide a guidance for a deeper theoretical 
understanding of hypernuclear dynamics and decay mechanisms. 

For a comprehensive review on the subject of these lectures we refer the reader to 
Ref. [5] and references therein. 



2. — Weak decay modes of A— hypernuclei 

In the production of hypernuclei, the populated state may be highly excited, above 
one or more threshold energies for particle decays. These states are unstable with respect 
to the emission of the hyperon, of photons and micleons. The spectroscopic studies of 
strong and electromagnetic de-excitations give information on the hypernuclear structure 
which are complementary to those we can extract from excitation functions and angular 
distributions studies. Once the hypcrmiclcus is stable with respect to electromagnetic 
and strong processes, it is in the ground state, with the hyperon in the Is level, and can 
only decay via a strangeness-changing weak interaction, through the disappearance of 
the hyperon. 
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with a Hfetime rf = ^rA*'*' = 2.632 x IQ-^" sec. 

Semi-leptonic and weak radiative A decay modes have negligible branching ratios: 

A ^ n-y (B.R. = 1.75 x IQ-^) 

pn-j (B.R. = 8.4 X IQ-^) 

pe-Ve (B.R. = 8.32 x IQ-^) 

PIJ,-V^ (B.R. = 1.57 X 10-*) 

and will not be considered here. 

The A hyperon is an isospin singlet (Ja = 0), while the ttA'' system can be either 
in / = 1/2 or in 7 = 3/2 isospin states. The customary angular momentum coupling 
implies: 
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Hence the ratio of amplitudes for A7 =1/2 transitions yields: 
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while a AI = 3/2 process should give: 
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Experimentally the above ratio turns out to be: 




Exp 



~ 1.78 



which is very close to 2 and strongly suggests the AI =1/2 rule on the isospin change. 
Prom the above considerations and from analyses of the A polarization observables it 
follows that the measured ratio between A7 = 1/2 and A7 = 3/2 transition amplitudes 
is very large: 



The AI =1/2 rule is based on experimental observations but its dynamical origin is not 
yet understood on theoretical grounds. It is also valid for the decay of the S hyperon 
and for pionic kaon decays (namely in non-leptonic strangeness-changing processes). 
Actually, this rule is slightly violated in the A free decay, and it is not clear whether it 
is a universal characteristic of all non -lcptonic processes with AS* 7^ 0. The A free decay 
in the Standard Model can occur through both AI — 1/2 and AI — 3/2 transitions, 
with comparable strengths: an s quark converts into a u quark through the exchange of 
a W boson. Moreover, the effective 4-quark weak interaction derived from the Standard 
Model including perturbative QCD corrections gives too small 1^1/2/^3/2! ratios (~ 3^4, 
as calculated at the hadronic scale of about 1 GeV by using renormalization group 
techniques [6]). Therefore, non-perturbative QCD effects at low energy (such as hadron 
structure and reaction mechanism), which are more difficult to handle, and/or final state 
interactions could be responsible for the enhancement of the AI =1/2 amplitude and/or 
the suppression of the A7 = 3/2 amplitude(^). 

The Q-valuc for frce-A mcsonic decay at rest is Q\ ~ ttia— mAr— ~ 40 McV. Then, 
taking into account energy-momentum conservation, niA — \/p^ + + ^/p^ + mj^ in 
the center-of-mass system and the momentum of the final nucleon turns out to be p ~ 100 
MeV. Inside a hypernucleus, the binding energies of the recoil nucleon (Bjv — —8 MeV) 
and of the A {Ba > — 27 MeV) tend to further decrease Qa [Qa, bound = Qa + Ba — Bjv] 
and hence p. 

As a consequence, in nuclei the A mesonic decay is disfavoured by the Pauli principle, 
particularly in heavy systems. It is strictly forbidden in normal infinite nuclear matter 



Al/2 

^3/2 



~ 30. 



(^) See, for example, a recent work based on the Instanton Liquid Model [7] 
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(where the Fermi momentum is kp ~ 270 MeV) , while in finite nuclei it can occur because 
of three important effects: 

1. In nuclei the hyperon has a momentum distribution (being confined in a limited 
spatial region) that allows larger momenta to be available to the final nucleon; 

2. The final pion feels an attraction by the medium such that for fixed momentum ^it 
has an energy smaller than the free one [u){q) < \/q'^ + m^], and consequently, due 
to energy conservation, the final nucleon again has more chance to come out above 
the Fermi surface. Indeed it has been shown [8, 9] that the pion distortion increases 
the mesonic width by one or two orders of magnitude for very heavy hypernuclei 
{A ~ 200) with respect to the value obtained without the medium distortion; 

3. At the nuclear surface the loc;al Fermi momentum is considerably smaller than k^, 
and the Pauli blocking is less effective in forbidding the decay. 

In any case the mesonic width rapidly decreases as the nuclear mass number A of the 
hypernucleus increases. 

The mesonic channel also gives information on the pion-nucleus optical potential 
since Fm = T^- + F^o is very sensitive to the pion self-energy in the medium: the 

latter is enhanced by the attractive P-wave 7r-nuclcus interaction and reduced by the 
repulsive iS'-wave one. Evidence for a central repulsion in the A-nucleus mean potential 
was obtained from the mesonic decays of s-shell hypernuclei [10, 11]. 

2'2. Non-mesonic decay. - In hypernuclei the weak decay can occur through processes 
which involve a weak interaction of the A with one or more nucleons. Sticking to the 
weak hadronic vertex A ttN, when the emitted pion is virtual, then it will be absorbed 
by the nuclear medium, resulting in a non-mesonic process of the following type: 

(1) An nn (F„) , 

(2) Ap ^ np (Fp) , 

(3) ANN nNN (F2) . 

The total weak decay rate of a A-hypernucleus is then: 

Tt = Tm + Tnm, 

where: 

Tm = r^- + r^o, Tnm = Ti + F2, Fi = F„ + Fp, 

and the lifetime is r = /i/Ft. The channel (3) can be interpreted by assuming that the 
pion is absorbed by a pair of nucleons, correlated by the strong interaction. Obviously, 
the non-mesonic processes can also be mediated by the exchange of more massive mesons 
than the pion (see figure 1). 
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Fig. 1. - One-nucleon (a) and two-nucleon (b) induced A decay in nuclei. 



The non-mesonic mode is only possible in nuclei and, nowadays, the systematic study 
of the hypemuclear decay is the only practical way to get information on the weak 
process AA^ NN (which provides the first extension of the weak AS* = NN NN 
interaction to strange baryons), especially on its parity-conserving part, which is masked 
by the strong interaction in the weak NN — > NN reaction. 

The final nucleons in the non-mesonic processes emerge with large momenta: dis- 
regarding the A and micleon binding energies and assuming the available energy Q = 
iTiA — ruN — 176 MeV to be equally splitted among the final nucleons, it turns out that 
Pn 420 MeV for the one-nucleon induced channels [Eqs. (1), (2)] and pn — 340 MeV 
in the case of the two-nucleon induced mechanism [Eq. (3)]. Therefore, the non mcsonic 
decay mode is not forbidden by the Pauli principle: on the contrary, the final nucleons 
have great probability to escape from the nucleus. The non-mesonic mechanism dom- 
inates over the mesonic mode for all but the s-shell hypernuclei. Only for very light 
systems the two decay modes are competitive. 

Since the non-mesonic channel is characterized by large momentum transfer, the 
details of the hypernuclear structure do not have a substantial influence (then providing 
useful information directly on the hadronic weak interaction). On the other hand, the 
NN and AA^ short range correlations turn out to be very important. 

It is interesting to observe that there is an anticorrelation between mesonic and non- 
mesonic decay modes such that the experimental lifetime is quite stable from light to 
heavy hypernuclei [12, 13], apart from some fluctuation in light systems because of shell 
structure effects: ta = (0.5 -i- 1)tjI'^'^. Since the mesonic width is less than 1% of the 
total width for A > 100, the above consideration implies that the non-mesonic rate is 
rather constant in the region of heavy hypernuclei. 

This can be simply understood from the following consideration. If one naively as- 
sumes a zero range approximation for the non-mesonic weak interaction AA" nN , then 
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Fig. 2. - Qualitative behaviour of masonic, non-mesonic and total decay widths as a function 
of the hypernuclear baryonic number A + 1. 



Ti is proportional to the overlap between the A wave function and the nuclear density: 



where the A wave function i/ja (nuclear density pa) is normalized to unity (to the nuclear 
mass number A). This overlap integral increases with the mass number and reaches 
a constant value: by using, e.g., A harmonic oscillator wave functions (with frequency 
oj adjusted to the experimental hyperon levels in hypernuclei) and Fermi distributions 
for the nuclear densities, we find ri(\2c)/ri(i"Spb) ~ 0.56, while Ti is 90 % of the 
saturation value for A ~ 65. In figure 2 the qualitative behaviour of mesonic, non- 
mesonic and total widths as a function of the nuclear mass number A is shown. 

For A < 11 the experimental data are quite well fitted by rNM/r^'''' — O.IA: Fi 
is proportional to the number of AN pairs. A, as it is expected from the above simple 
description, where we neglect the contribution of F2. However, Fi saturates when the 
radius of the hypernucleus becomes sensitively larger than the range of the AA^ nN 
interaction. For a more quantitative explanation it will be important to collect data with 
good precision. Yet, from the available data one can roughly say that the long distance 
component of the A A'' nN transition has a range of about 1.5 fm and corresponds, as 
we expect, to the one-pion-exchange component of the interaction. 
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2'3. The Tn/Tp puzzle. Nowadays, an important question concerning the weak decay 
rates is the longstanding disagreement between theoretical estimates and experimental 
determinations of the ratio Tn/Tp between the neutron- and the proton-induced decay 
widths. 

This problem will be extensively discussed in Section 5. However, it is worth re- 
calling here that, up to short time ago, all theoretical calculations appeared to strongly 
underestimate the available central data measured in several hypernuclei: 

Until recently the data were quite limited and not precise because of the difficulty in 
detecting the products of the non-mesonic decays, especially the neutrons. Moreover, 
the experimental energy resolution for the detection of the outgoing miclcons does not 
allow to identify the final state of the residual nuclei in the processes -^"^Z + nn 

and ^~^(Z — 1) + np. As a consequence, the measurements supply decay rates 

averaged over several nuclear final states. 

In the one-pion-exchangc (OPE) approximation, by assuming the A/ ~ 1/2 rule in 
the A TT~p and A tt^ti free couplings, the calculations (which will be reported 
later) give small ratios, in the range 0.05 -t- 0.20 for all the considered systems. However, 
as we shall see in Section 4, the OPE model with A/ =1/2 couplings has been able to 
reproduce the one-body stimulated non-mesonic rates Fi = F^+Fp for light and medium 
hypernuclei. Hence, the problem seems to consist in overestimating the proton-induced 
and underestimating the neutron induced transition rates. 

In order to solve this puzzle (namely to explain both F„ + Fp and F„/Fp), many 
attempts have been made up to now, mainly without success. We recall the inclusion in 
the AN — > nN transition potential of mesons heavier than the pion (also including the 
exchange of correlated or uncorrelated two pions) [14-18], the inclusion of interaction 
terms that explicitly violate the AI = 1/2 rule [19] and the description of the short 
range baryon-baryon interaction in terms of quark degrees of freedom [20,21], which 
automatically introduces A7 = 3/2 contributions. 

3. — Theoretical models for the decay rates 

We illustrate here the theoretical approaches which have been utilized for the formal 
derivation of A decay rates in nuclei. We discuss first the general features of the approach 
used for direct finite nucleus calculations. It is usually called Wave Function Method 
(WFM), since it makes use of shell model nuclear and hypernuclear wave functions (both 
at hadronic and quark level) as well as pion wave functions generated by pion-nucleus 
optical potentials. Then we consider the Polarization Propagator Method (PPM), which 
relies on a many-body description of the hyperon self-energy in nuclear matter. The 
Local Density Approximation allows then one to implement the calculation in finite 
nuclei. Finally, a microscopic approach, based again on the PPM, is shortly sketched: in 
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this case the full A self-energy is evaluated on the basis of Feynman diagrams, which are 
derived, within a functional integral approach, in the framework of the so-called bosonic 
loop expansion. 

3'1. Wave Function Method: mesonic width. - The weak effective Hamiltonian for 
the A — > nN decay can be parameterized in the form: 

(4) HZn = iGmli^^iA + B^^)t ■ ^V-a, 

where the values of the weak coupling constants G — 2.211 x 10^^/mJ, A — 1.06 and 
B = —7.10 are fixed on the free A decay. The constants A and B determine the strengths 
of the parity violating and parity conserving A nN amplitudes, respectively. In order 
to enforce the A/ = 1/2 rule (which fixes r^°7rJ^^o°° = 2), in Eq. (4) the hyperon is 
assumed to be an isospin spurion with I = 1/2, 1^ = —1/2. 

In the non-relativistic approximation, the free A decay width r5('°° = T^^'i'^ + F^o'' is 
given by: 



Fl- = c^iGmir I ^^Jl^^^ 2.5[m^ - - E^] [s^ + 



2 I _ ;?2 

ni: 



where c„ = 1 for F^o and = 2 for F^- (expressing the A/ = 1/2 rule), S = A, 
P = mTrB/{2mN), whereas En and u){q) are the total energies of nucleon and pion, 
respectively. One finds the well known result: 

F^- = c„(Gm2)2 J-™£:^iL ( + flq^ 



27r rriA V ml^^-'^- 

which reproduces the observed rates. In the previous equation, qc.m. — 100 MeV is the 
pion momentum in the center -of- mass frame. 

In a finite nucleus approach, the mesonic width Fm = F^- + F^o can be calculated 
by means of the following formula [8, 9]: 



= c^iGmir J 2.6[E. - .(g) - E.] 



df(j)A{f)(l)„{q,r)(l)%{r) 



2 p2 



2^ 

df(l)A{r)V(pn{q, r)4'*N{r) 



where the sum inns over non-occupied nuclconic states and i?A is the hyperon total 
energy. The A and nucleon wave functions (/)a and (pN are obtainable within a shell 
model. The pion wave function corresponds to an outgoing wave, solution of the 
Klein-Gordon equation with the appropriate pion-nucleus optical potential 14pt: 



-ml- 2u;l/opt(r) + Vc{f)f] 0.(«, f) 
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where Vc{r) is the nuclear Coulomb potential and the energy eigenvalue ui depends on 
Q- 

Different calculations [8, 9] have shown how strongly the mesonic decay is sensitive to 
the pion-nucleus optical potential, which can be parameterized in terms of the nuclear 
density, as discussed in Refs. [9], or evaluated microscopically, as in Ref. [8]. 

3"2. Wave Function Method: non-mesonic width. - Within the one-meson-exchange 
(OME) mechanism, the weak transition AA'' — > nN is assumed to proceed via the media- 
tion of virtual mesons of the pseudoscalar (tt, t] and K) and vector (p, co and K*) octets 

[14,15] (sec Fig 1). 

The fundamental ingredients for the calculation of the AN nN transition within 
a OME model are the weak and strong hadronic vertices. The AwN weak Hamiltonian 
is given in Eq. (4). For the strong NNn Hamiltonian one has the usual pseudoscalar 
coupling: 

gNN-K being the strong coupling constant for the NN-k vertex. In momentum space, the 
non-relativistic transition potential in the OPE approximation is then: 

.r n 2 9NN-n f . B ^ \ (72 ■ q ^ ^ 

2mjv V 2m J q-^ + 

where m = {m\ + mN)/'2 and q is the momentum of the exchanged pion. 

Due to the large momenta (~ 420 McV) exchanged in the AA'' nN transition, 
the OPE mechanism describes the long range part of the interaction, and more massive 
mesons are expected to contribute at shorter distances. 

Non-trivial difficulties arise with the heavier mesons, since their weak couplings in the 
AA^ vertex are not known experimentally. For example, if one includes in the calculation 
the contribution of the p-meson, the weak ANp and strong NNp Hamiltonians give rise 
to the following p-meson transition potential: 



V,{q) = Gml 



9nnp<^ 7 —{<^i xq)-{a2Xq) 



, ■^idNNp'^ 9NNp) ,-, ^ ^ - 
+1 ^ '-((Ti X (72) • 5 

2mm 



Tl ■ T2 

q'^ + m?p 



where the weak coupling constants a, /3 and e must be evaluated theoretically and turn 
out to be quite model-dependent. 

The most general OME potential accounting for the exchange of pseudoscalar and 
vector mesons can be expressed through the following decomposition: 

(5) V{r) = Y.Vm{7^ = 5]5]C(r)6«(rV™, 
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where m = n, p, K, K*, to, rj; the spin operators O" are (PV stands for parity-violating): 



1 



0"(r) = < 



(Tl • (72 



(72 • r 

(iTi X CT2) • f 



central spin-independent, 
central spin-dependent , 



•512 (r) = 3{ai ■ r){ff2 ■ r) — a\ ■ $2 tensor, 



PV for pseudoscalar mesons, 
PV for vector mesons, 



whereas the isospin operators are: 



isoscalaxs mesons (ry, w), 
isovector mesons (tt, p), 



Jinear combination of 1 and ti • T2 isodoublet mesons (if, K*). 

For details concerning the potential (5), see Ref. [15, 17]. 

Assuming the initial hypcrnucleus to be at rest, the one-body induced non-mesonic 
decay rate can then be written as: 



(6) 



(27r)3J (27r)3 

where 5(E.C.) stands for the energy conserving delta function: 



Moreover: 



M{puP2) = r\N{px)N{p2)\Tan^nn\^ h) 

is the amplitude for the transition of the initial hypernuclear state \^h) of mass mn into 
a final state composed by a residual nucleus \^ r) with energy Er and an antisymmetrized 
two nucleon state \N{pi)N{'p2)), pi and p2 being the nuclcon momenta. The sum ^ in 
Eq. (6) indicates an average over the third component of the hypernuclear total spin and 
a sum over the quantum numbers of the residual system and over the spin and isospin 
third components of the outgoing miclcons. Customarily, in shell model calculations 
the weak-coupling scheme is used to describe the hypernuclear wave function \'^h), the 
nuclear core wave function being obtained through the technique of fractional parentage 
coefficients [15]. The many-body transition amplitude A^(pi,P2) is then expressed in 
terms of two- body amplitudes {NN\V\AN) of the OME potential of Eq. (5). 
Two merits of the WFM must be remarked: 
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• Since the A decays from an orbital angular momentum / = state, in the non- 
mesonic decay rate one can easily isolate the contributions of neutron- and proton- 
induced transitions [15], and the Tn/Tp ratio can be directly evaluated. 

• The nN final state interactions and the AA'' correlations (which are absent in 
an independent particle shell model) can also be implemented in the calculation 
[15,17]. 

3 '3. Polarization Propagator Method and Local Density Approximation. - The hyper- 
nuclear decay rates can be studied by using the Polarization Propagator Method [22] to 
evaluate the A self-energy inside the nuclear medium. The polarization propagator is 
conveniently calculated for a homogeneous system (nuclear matter) , within the Random 
Phase Approximation (RPA) and eventually accounting for additional correlations. The 
calculation can then be extended to finite nuclei via the Local Density Approximation 
(LDA). 

This many body technique provides a unified picture of the different decay channels 
and it is equivalent to the WFM [23] (in the sense that it is a semiclassical approxima- 
tion of the exact quantum mechanical problem). Obviously, for the mesonic rates the 
WFM is more reliable than the PPM in LDA, since Fm is rather sensitive to the shell 
structure of the hypernucleus. On the other hand, the propagator method in LDA offers 
the possibility of calculating the hypernuclear decay rates over a broad range of mass 
numbers, while the WFM is hardly exploitable for medium and heavy hypernuclei. 

To calculate the hypernuclear width one needs the imaginary part of the A self-energy: 

(7) rA = -2ImSA, 
which, in the non-relativistic limit, reads: 

(8) SA(fc) = HGmlfJ ^ (^S^ + ^q^^ Fl{q)GN{k - q)G.{q). 
The nucleon and pion propagators in nuclear matter are, respectively: 

(9) G.(p) = t(l^^|-t^ . + I ^^i) 



Po - En{p) -VN + ie Po- En(p) -Vn -ie' 
and: 

(10) G^{q) 



ql-q^ -ml-Y.l{qy 



The above form of the non-relativistic nucleon propagator refers to a non-interacting 
Fermi system but includes corrections due to Pauli principle and an average binding. 

In the previous equations, p = {p(),p) and q = {qo,q} denote four vectors, kp is the 
Fermi momentum, E^^ is the nucleon total free energy, Vn is the nucleon binding energy 




(f) (g) (h) 

Fig. 3. - Lowest order terms for the A self-energy in nuclear matter. The meaning of the various 
diagrams is explained in the text. 



(which is density-dependent), and S* is the pion proper self-energy in nuclear matter. 
This quantity has been carefully evaluated within several many body frameworks and 
includes the (strong) coupling of the pion to particle-hole {p~h) states, collective RPA 
states and more complicated nuclear correlated states (e.g., the two-particle two-hole 
states). A monopole form factor F„{q) describing the hadronic structure of the ttAA'' 
vertex is also included in Eq. (8). 

We note here that the parity-conserving term {1 = 1 term) in Eq. (8) contributes 
only about 12% of the total free decay width. However, the P-wave interaction becomes 
dominant in the nuclear non-mesonic decay, because of the larger exchanged momenta. 

In Fig. 3 wc show the lowest order Feynman diagrams for the A self-energy in nuclear 
matter. Diagram (a) represents the bare self-energy term, including the effects of the 
Pauli principle and of binding on the intermediate nucleon. In (b) and (c) the pion 
couples to a p h and a A i pair, respectively. Diagram (d) is an insertion of S-wave 
pion self-energy at lowest order. In diagram (e) we show a 2p-2h excitation coupled to 
the pion through S-wa.ve ttN interactions. Other 2p-2h excitations, coupled in P-wave, 
are shown in (f) and (g), while (h) is a RPA iteration of diagram (b). 

In Eq. (8) there are two different sources of imaginary part. The analytical structure 
of the integrand allows the integration over qq [24]. After performing this integration, 
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an imaginary part is obtained from the (renormalized) pion-nucleon pole and physically 
corresponds to the mesonic decay of the hyperon. Moreover, the pion proper self-energy 
E* (g) has an imaginary part itself for {qq, q) values which correspond to the excitation of 
p-h, A-h, 2p-2h, etc states on the mass shell. By expanding the pion propagator 0,^(0) 
as in Fig. 3 and integrating Eq. (8) over qo, the nuclear matter A decay width of Eq. (7) 
becomes [24]: 

TA{k, p) = -6{Gmlf J k-q\ -kF)e{ko - E^ik - q) - Vn) 

(11) H<l)]<io=ko-E^{k-g)-V^ ' 



where: 



(12) 



, Pliq)UL{q) P^{q)UT{q) 
^1 - VL{q)UL{q) 1 - VT{q)UT{q) ' 



In Eq. (11) the first 6 function forbids intermediate micleon momenta smaller than the 
Fermi momentum, while the second one requires the pion energy q^ to be positive. More- 
over, the A energy, fco = Ef^(k) + Va, contains a phenomenological binding term. With 
the exception of diagram (a), the pion lines of Fig. 3 have been replaced, in Eq. (12), 
by the effective interactions S, Pl, Pt,Vl, Vt {L and T stand for spin-longitudinal and 
spin-transverse, respectively), which include tt- and p-exchange modulated by the effect 
of short range repulsive correlations. The potentials Vl and Vr represent the (strong) p-h 
interaction and include a Landau parameter 9', which accounts for the short range repul- 
sion, while S, Pl and Pt correspond to the lines connecting weak and strong hadronic 
vertices and contain another Landau parameter, g'j^, which is related to the strong A A'' 
short range correlations. 
Furthermore, in Eq. (12): 

Gl{q) = — ^2 2"' 

q^-q^ - mi 

is the free pion propagator, while UL{q) and UT{q) contain the Lindhard functions for 
p-h and A-Ji excitations [25] and also account for the irreducible 2p-2h polarization 
propagator: 

(13) Ul,Aq) = U^\<l) + U'^\<l) + Uf^ri^l)- 
They appear in Eq. (12) within the standard RPA expression: 

rr^^T^rr TT " ^L{T)iq) + U L(T){q)VL{T){q)U L{T){q) + 

1 - Vz,(T)(?)J7i(T)(«) 

(14) +UL{T){q)VL{T){q)UL(T){q)VL{T){q)UL{T){q) + ■■■ 
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The decay width (11) depends both exphcitly and through Ul,t{q) on the nuclear matter 
density p = 2fc|./37r^. 

The Lindhard function JJp'' {U^^) is given by: 

(15) UP^^'^^\q) = -M I ^G%{p)G%^^^{p + q), 

with the free nucleon and Delta propagators: 

9{\p\-kF) , 0{kF-\p\) 



(16) G%{p) 



Po -Tn{p}+ ie po -Tn{p) -ie' 



(17) Glip) = ^ ^ 



9 - 7a(p) - 6M^N + i rA/2 ' 



where ^^(A) is the nucleon (Delta) kinetic energy, Fa the free A width and (5Majv = 
TOA — TO AT. We remind that UP^^iq) and W^'^iq) can be analytically evaluated in nuclear 
matter. 

For the evaluation of U'^'^ wc discuss two different approaches. In Refs. [26, 27] 
a phenomenological parameterization was adopted: this consists in relating to 
the available phase space for on-shell 2p-2h excitations in order to extrapolate for off- 
mass shell pions the experimental data of P~wave absorption of real pions in pionic 
atoms. In an alternative approach [28], U"^^ is evaluated microscopically, starting from 
a classification of the relevant Feynman diagrams according to the so-called bosonic loop 
expansion, which is obtained within a functional approach. 

Wc notice here, in connection with Eqs. (8) and (10), that the full (proper) pion 
self-energy 

s;(g) = s(.^)*(g) + E(^)*(g), 

contains a P-wave term, which is related to the spin-longitudinal polarization propagator 
UL{q) according to: 



(18) 4'')*(g) = 



f^Fl{q)UL{q) 
l-AgL{q)UL{q) 



the Landau function (q) being given in the appendix of Rcf. [5] . E* (q) also contains 
an S-wave term, which, by using the parameterization of Ref. [29], can be written as: 

(19) Ei^^*{q)=-in(l + ^]bop. 
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Fig. 4. - Schematic representation of the total A self-energy in terms of the 7r-self energy and 
polarization insertions; in the vertexes w = weak, s = strong. 



The parameter bo is usually taken from the phenomenology of pionic atoms (see, for 
example, Ref. [30]). The function T,^^* is real (constant and positive), therefore it 
contributes only to the mesonic decay [diagram (d) in Fig. 3 is the relative lowest order]. 
On the contrary, the P-wave self-energy is complex and attractive: Re eI^'' *{q) < 0. It 
contributes to all the decay channels, in particular to the two-body induced non-mesonic 
decay width. A schematic picture of the relation between the A and tt self-energies 
is illustrated in Fig. 4, where the polarization propagator insertion on the pionic line 
summarizes all contributions explicitly shown in Fig. 3 and many others, up to infinite 
order. 

The propagator method provides a unified picture of the decay widths. A non- 
vanishing imaginary part in a self-energy diagram requires placing simultaneously on- 
shell the particles of the considered intermediate state. For instance, diagram (b) in 
Fig. 3 has two sources of imaginary part. One comes from cut 1, where the nucleon and 
the pion are placed on shell. This term contributes to the mesonic channel: the final 
pion eventually interacts with the medium through a p-h excitation and then escapes 
from the nucleus. Diagram (b) and further iterations lead to a renormalization of the 
pion in the medium which may increase the mesonic rate even by one or two orders of 
magnitude in heavy nuclei [8,9,24]. The cut 2 in Fig. 3(b) places a nucleon and a p-h 
pair on shell, so it is the lowest order contribution to the physical process AA^ nN. 
Analogous considerations apply to all the considered diagrams. 

A couple of remarks may be useful, here, for those who are not familiar with the 
language of many-body theory. 

1. The one body induced process AN nN (for example via the exchange of a pion) 
translates, inside the nuclear medium, into the creation of a particle-hole pair, the 
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(particle) 
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Fig. 5. - One-body induced process in the language of particle-hole states. 



hole representing the state of the initial nucleon, which turns out to be vacant, in 
the medium, after the weak interaction with the A. This is schematically illustrated 
in Fig. 5. 

2. The relation between the one-body induced decay width, Fi, and the imaginary 
part of the A self-energy, say, of diagram (b) in Fig. 3, is represented in the upper 
panel of Fig. 6, where the modulus square of the decay amplitude is pictorially 
related to the Imaginary part of the considered diagram, with a cut putting on- 
shell the 2 nucleon, 1 hole intermediate state. Analogously, the two-body induced 
decay width is related to the Imaginary part of diagrams like, e.g., Fig. 3(f), the 
cut being placed on the 2p 2h intermediate state of the polarization propagator. 
This is represented in the lower panel of Fig. 6. 

In order to evaluate the various contributions to the width stemming from Eq. (11), 
it is convenient to consider all the intervening free meson propagators as real. Then the 
imaginary part of (12) will develop the following contributions: 

™ 1 - VL,T{q)UL,Tiq) I 1 - Vl,t{q)Ul,t{q) P 

The three terms in the numerator of Eq. (20) can be interpreted as different decay 
mechanisms of the hypcrnuclcus. The term proportional to ImU^^ provides the one— 
nucleon induced non-mesonic rate, F,. There is no overlap between Im (g) and the 
pole Qo = u}{q) in the (dressed) pion propagator G.„{q): thus the separation of the mesonic 
and one-body stimulated non-mesonic channels is unambiguous. 

Further, Im U^^ accounts for the A ttN decay width, thus representing a (small) 
contribution to the mesonic decay. 

The third contribution of Eq. (20), proportional to ImU^^^, intervenes in a wide 
kinematic range, in which the above mentioned cuts put on the mass shell not only the 
2p-2h lines, but possibly also the pionic lines. Indeed, the renormalized pion pole in 
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Im 



2p-2h polar, prop. 




( strong) ph-interaction 



Fig. 6. - Schematic relation between the one-body (upper panel) and two-body (lower panel) A 
decay amplitudes and the Imaginary part of specific contributions to the A self-energy. 



Eq. (10) is given by the dispersion relation: 



m2-ReE*[a;(g-),g]=0, 



with the constraint: 



u!{q) = ko- EN{k -q)-VN. 



At the pion pole, \m.U'^^ ^ 0, thus the two body induced non-mesonic width, 
cannot be disentangled from the mesonic width, Fm- In other words, part of the decay 
rate calculated from Im is due to the excitations of the renormalized pion and gives 
in fact Fm, with the exception of the mesonic contribution originating from IvnU^^, 
which is, however, only a small fraction of Fm- 

In order to separate Fm from F2, in the numerical calculation it is convenient to eval- 
uate the mesonic width by adopting the following prescription. We start from Eq. (11), 
setting: 



(21) 



a{q)=aM{q) = \ S"^ + 



P2 



mt 



F'{q)GAq), 



and omitting Im S* in Gtt (which corresponds to assume Im U^'^ = Im t/ 



Ah 



ImC/^f^'' 
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0). Then Imauiq) '^^^Y accounts for the (real) contribution of the pion pole: 
Im G^{q) = -ttS [q^ - f -ml -Re S;(g)] . 

Once the mesonic decay rate is known, one can calculate the three-body non-mesonic 
rate by subtracting Fm and Fi from the total rate Ft, which one gets via the full 
expression for a (g) [Eq. (12)]. 

Using the PPM, the decay widths of finite nuclei can be obtained from the ones 
evaluated in nuclear matter via the LDA: in this approach the Fermi momentum is made 
r-dependent (namely a local Fermi sea of nuclcons is introduced) and related to the 
nuclear density by the same relation which holds in nuclear matter: 



(22) 



1/3 



More specifically, one usually assumes the nuclear density to be a Fermi distribution: 



PA{r) 



1 + exp 



r - R{A) 



with R{A) = 1.12A1/3-0.86A-1/3 fm, a = 0.52 fm and po = a\ ^irR^iA) ( 1 



R(A) 

Moreover, the nucleon binding potential Vn also becomes i — dependent in LDA. In 
Thomas-Fermi approximation one assumes: 



eF(r) + Viv(r) 



2mjv 



+ Fjv(fO = 0. 



With these prescriptions one can then evaluate the decay width in finite nuclei by 
using the semiclassical approximation, through the relation: 



(23) 



rA(fc) = J 



df|VA(r)l^rA k,p{f) 



where ipA is the appropriate A wave function and TA[k, p{f)] is given by Eqs. (11), (12) 
with a position dependent nuclear density. This decay rate can be regarded as the k- 
component of the A decay rate in the nucleus with density p{r). It can be used to 
estimate the decay rates by averaging over the A momentum distribution ItpAik)]"^- One 
then obtains the following total width: 



(24) 



/ 



dk\Mk)\''TA{k), 



which can be compared with the experimental results. 
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3'3.1. Phenomenological 2p-2h propagator. Coming to the phenomenological evalu- 
ation of the 2p-2h contributions in the A self-energy, the momentum dependence of the 
imaginary part of U'fj'^'^ can be obtained from the available phase space, through the 
following equation [26]: 

(25) ImC/|P^''(go,g;p) = ^ Im ^^j,'' (m^ , 0; p^s ) , 

P(m^,0;peflF) 

where peS = 0.75/9. By neglecting the energy and momentum dependence of the p- 

h interaction, the phase space available for on- shcU 2p-2h excitations [calculated, for 
simplicity, from diagram 3(e)] at energy-momentum {qo,q) and density p turns out to 
be: 

P{qo, q; p) (X 1 ^ Im U^^ (| + k; p) Im C/f'^ (| - k; p) 

X.(|..o).(|-.o). 

In the region of {qo,q) where the p h and A h excitations are off- shell, the relation 
between u'^p'^^ and the P wave pion-nucleus optical potential Vopt is given by [see also 
Eq. (18), in the language of pion self-energy]: 

(26) = 2qoVoM. 

1 - ^9L{q)UL{q) 
mi 

At the pion threshold Vopt is usually parameterized as: 

(27) 2qoVopt{qo ^m^,qc^ 0; p) = -Anq'^p'^Co, 

where Co is a complex number which can be extracted from experimental data on pionic 

atoms. By combining Eqs. (26) and (27) it is possible to parameterize the proper 2p-2h 
excitations in the spin-longitudinal channel through Eq. (25), by setting: 

(28) ^^Ilp^iqo - m,, qc^ 0)C/f '''(go - m^, g ~ 0; p) = -Anq^p^C^. 

The relation between Co and Cq is fixed on the basis of the same RPA expression for 
the polarization propagator contained in Eq. (26); hence the value of Cq also depends 
on the correlation function g^- From the analysis of pionic atoms data made in Ref. [31] 
and taking g' = gL{0) = 0.615, one obtains: 



Co = (0.105 + z0.096)/m^ 
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The spin transverse component of U^p'^'^ is assumed to be equal to the spin-longitudinal 
one, U^^^'^ — [/^^^'*, and the real parts of JJ^^'^ and U^^"^^ are considered constant [by 
using Eq. (28)] because they arc not expected to be too sensitive to variations of qo and 
q. The assumption J/^^^'' = jj'^p'^''- jg not a priori a good approximation, but it is the 
only one which can be employed in the present phenomenological description. Yet, the 
differences between [/^p^'* and U^'^^ can only mildly change the partial decay widths [5]: 
arc summed to f^^, which gives the dominant contribution. Moreover, for 
jj2p2h _ jj^p'^h ^Y^^ transverse contribution to r2 [fourth term in the right hand side of 
Eq. (12)] is only about 16% of r2 (namely 2 3% of the total width) in medium-heavy 
hypernuclei. 

The simplified form of the phenomenological 2p-2h propagator, together with the 
availability of analytical expressions for U^'^ and U^~'^, makes this approach particularly 
suitable for employing the above mentioned LDA. 

3'3.2. Functional approach to the A self-energy. In alternative to the above men- 
tioned phenomenological approach for the two-body induced decay width, we briefly 
discuss here a truly microscopic approach. Indeed the most relevant Fcynman diagrams 
for the calculation of the A self-energy can be obtained in the framework of a functional 
method: following Ref. [28], one can derive a classification of the diagrams according to 
the prescription of the so-called bosonic loop expansion. 

The functional techniques can provide a theorc;tically founded derivation of new classes 
of Feynmann diagram expansions in terms of powers of suitably chosen parameters. For 
example, the already mentioned ring approximation (a subclass of RPA) automatically 
appears in this framework at the mean field level. This method has been extensively ap- 
plied to the analysis of different processes in nuclear physics [32 34] . Here it is employed 
for the calculation of the A self-energy in nuclear matter, which can be expressed through 
the nuclear responses to pseudoscalar-isovector and vector-isovector fields. The polar- 
ization propagators obtained in this framework include ring dressed meson propagators 
and almost the whole spectrum of 2p-2h excitations (expressed in terms of a one-loop 
expansion with respect to the ring-dressed meson propagators), which are required for 
the evaluation of 

Let us first consider the polarization propagator in the pionic (spin-longitudinal) 
channel. To illustrate the procedure, it is useful to start from a Lagrangian describing a 
system of nucleons interacting with pions through a pseudoscalar-isovector coupling: 

(29) jC^n = - mjv)V + ^^^^ • 9^"$ - ^mlP - ii>f^l) ■ $, 

where V {4') is the nucleonic (pionic) field and F = g^^r {g = 2fT^mN /m.^) is the spin- 
isospin matrix in the spin-longitudinal isovector channel. We remind the reader that in 
the calculation of the hypernuclear decay rates one also needs the polarization propagator 

in the transverse channel [see Eqs. (11) and (12)]: hence, one has to include in the model 
other mesonic degrees of freedom, like the p meson. Since the bosonic loop expansion is 
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characterized by the topology of the diagrams, this is relatively straightforward and the 
same scheme can be easily applied to mesonic fields other than the pionic one. 

The generating functional, expressed in terms of Feynman path integrals, associated 
with the Lagrangian (29) has the form: 



(30) 



V 



'ip,ip,(l) 



exp 



dx 



C^n{x) - iii){x)f^{x) ■ 0{x) 



where a classical external field ip with the quantum numbers of the pion has been in- 
troduced (here and in the following the coordinate integrals are 4-dimensional). All the 
fields in the functional integrals have to be considered as classical variables, but with the 
correct commuting properties (hence the fermionic fields are Grassman variables). 

The physical quantities of interest for the problem are then deduced from the gener- 
ating functional by means of functional differentiations. In particular, by introducing a 
new functional Zn such that: 



(31) 



Z[0] = exp{iZc[0]}, 



the spin longitudinal, isovcctor polarization propagator turns out to be the second func- 
tional derivative of Zc with respect to the source (p of the pionic field: 



(32) 



S^Z,[0] 
S(pi{x)S(pj{y) 



We note that the use of Zc instead of Z in Eq. (32) amounts to cancel the disconnected 

diagrams of the corresponding perturbative expansion (linked cluster theorem). From 
the generating functional Z one can obtain different approximation schemes according 
to the order in which the functional integrations are performed. 

For the present purpose, it is convenient to integrate Eq. (30) over the nucleonic 
degrees of freedom first. Introducing the change of variable cf) ^ (j) — one obtains: 



(33) 



Z[(p] = exp 



1/ 



dx dy ip{x) • G° \x - y)(p{y) 



X / V 



exp<^z / dxdy \ij{x)Gj/{x - y)ip{y) 



I 

l${x)-Gl-\x-y) {${y) + 20{y)) 



where the integral over [tp,^^] is gaussian: 



exp<ji / dxdy'il;{x)G^^{x - y)'ip{y) \ = {det Gn) ■ 
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Hence, after multiplying Eq. (33) by the unessential factor detG^^ {G^ being the free 
nucleon propagator), which only redefines the normalization constant of the generating 
functional, and using the property detX = exp{TrlnX}, one obtains: 



(34) 



Z[0] = exp I ^ / dx dy (p{x) • G° \x - y)0{y) 



X / V 



exp 



In this expression we have introduced the bosonic effective action: 



(35) 

where(^): 
(36) 14 



dx dy \ -$ix) ■ Gl" (x - y) ${y) + 20{y) 



n=l 

i^Tr(rir,) / dxdy^''{x,y)4>i{x)c^M 
^ E (^^-'"■5^'^) / dxdydzlf{x,y,z)(p,{ 



x)cPj{y)Mz)+Oi,^^). 



and: 

(37) 
(38) 



-in°(x, y) = tG%{x - y)iG%{y - x), 
-ill°{x, y, z) = iG%{x - y)iG%{y - z)iG%(z - x), etc. 



The bosonic effective action (35) contains a term for the free pion field and also a highly 
non-local pion self-interaction l^-. This effective interaction is given by the sum of all 
diagrams containing one closed fermion loop and an arbitrary number of pionic legs. 
We note that the function in Eq. (37) is the free particle-hole polarization propagator. 
Moreover, the functions lf{x,y,...,z) are symmetric for cyclic permutations of the 
arguments. 

The next step is the evaluation of the functional integral over the bosonic degrees of 
freedom in Eq. (34). A perturbative approach to the bosonic effective action (35) does 

(^) Eq. (36) is a compact writing: for example, the n = 2 term must be interpreted as: 



^Ti (if ■ $G%y = ^ I dxdyTrJ2i^iG°Nix - y)irjG%iy - x)(t>iix)<l>j{y), 
where the trace in the right hand side acts on the vertices f, and so on. 
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not seem to provide any valuable results within the capabilities of the present computing 
tools and in Ref. [28] another approximation scheme, the semiclassical method, was 
followed. 

The lowest order of the semiclassical expansion is the stationary phase approximation 
(also called saddle point approximation in the Euclidean space): the bosonic effective 
action is required to be stationary with respect to arbitrary variations of the fields (^ij: 



"'-'eff 



l5(/)i(x) 



= 0. 



From the partial derivative of Eq. (35) one obtains the following equation of motion for 
the classical field 6: 



(39) 



{n+ml)cPi{x)= j dyGl'\x-y)^,{y) + 







5cj>i 


X) 



whose solutions are functional of the external source 0. 

It can be shown that in the saddle point approximation the polarization propagator 
coincides with the well known ring expression: 



n°(a;,y)+Tr(r2) j dudv\i^{x,u)G'^''^{u - v)n^{v,y) 



5ijW'^^{x,y), 



where: 



1 - TV (r|) GO no 



is the ring-dressed pion propagator. Hence, the ring approximation corresponds to the 
mean field level of the present effective theory. 

In the next step of the semiclassical method, the bosonic effective action is expanded 
around the mean field solution: 



S. 



eff 



1 



S.2 oB 



5(j)r{x)6(t)j{y) 



X [Ux) - <tP,{x)] [cj>^{y)-cjP.{y)]. 

Then, after performing the gaussian integration over ^, the generating functional (34) 
[up to an uninteresting normalization factor] reads: 



(40) 



Z [(p] = exp <^ iS^s 



- -Trln 
2 



1(2 oB 



5(i)r{x)5(j).j{y) 
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(a) 



(b) 
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(d) 





(e) 



(0 



Fig. 7. - Feynman diagrams for the polarization propagator of Eq. (41): (a) particle-hole; 
(b) exchange; (c) and (d) self-energy-type; (e) and (f) correlation diagrams. Only the first 
contribution to the ring expansion has been drawn. The dashed lines represent ring-dressed 
pion propagators. 



Lengthy but straightforward calculations [28] lead then to the following total polarization 
propagator: 



(41) 
where: 



n,,- {x, y) = Sij [U"^^{x, y) + nOBi^(ar, y)] 



Il°'''^ix,y) = J2T^iW I dudvG^^iu- v)Il\x,u,y,v) 

kl '' 

+ ^ TV (r^^r,) / du dv G™8(u _ y) [11° {x, u, V, y) + n\x, y, v, u)] 

kl 

+ j dudv dw ds Gf'^^iu - w)G';''^{v - s)T^\x, u, v) 
(42) X [Tr(r,r,r )n°(2/,s,«;)] . 



klmn 



The Feynman diagrams corresponding to Eqs. (41) and (42) are depicted in Fig. 7. 
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Diagram (a) represents the Lindhard function n°(x,y), which is just the first term of 
H™s{x,y). In (b) we have an exchange diagram (the thick dashed hnes representing 
ring-dressed pion propagators); (c) and (d) are self-energy diagrams, while in (e) and (f) 
we show the correlation diagrams of the present approach. The approximation scheme 
discussed here is also referred to as bosonic loop expansion (BLE). The practical rule 
to classify the Feynman diagrams according to their order in the BLE is to reduce to a 
point all its fcrmionic lines and to count the number of bosonic loops left out. In this 
case the diagrams (b)-(f) of Fig. 7, which correspond to n°^^ Eq. (42), reduce to a 
one-boson-loop (OBL). 

The polarization propagator of Eq. (42) is the central result of the present microscopic 
approach and it has been used [28] for the calculation of the A decay width in nuclear 
matter. Note that the model can easily include the excitation of baryonic resonances, 
by replacing the fermionic field with multiplets. The topology of the diagrams remains 
the same as in Fig. 7 but, introducing for example the A resonance, each fcrmionic line 
represents either a nucleon or a A, taking care of isospin conservation. Obviously this 
procedure substantially increases the number of diagrams. 

Moreover, being the BLE characterized by the topology of the diagrams, additional 
mesonic degrees of freedom together with phenomenological short range correlations can 
be included by changing the definition of the vertices in Eq. (42). In particular, the 
formalism can be applied to evaluate the functions Ul,t of Eq. (13), which are required 
in Eqs. (11), (12). In the OBL approximation of Eq. (41) and Fig. 7 one has to replace 
Eq. (12) with: 

(43) ai,) = (S^ + ^gA FUg)Gl{,) ' ^'^'^^^^^'^^ 



l~VL{q)Ui{q) 



l-VL{q)Ui{q) l-VT{q)Ui{q) 



+ S\q)+Pl{q) C/r^9) + 2P|(«)C/OB^(g), 



where 



while are evaluated from the diagrams 7(b)-7(f) using the standard Feynman 

rules. The normalization of these functions is such that UP^{x,y) = 4lP{x,y), II" being 
given by Eq. (37). One relevant difference between the OBL formula (43) and the RPA 
expression of Eq. (12) lies in the fact that in the former, to be consistent with Eq. (42), 
the 2p-2h diagrams (which contribute to U^^) are not RPA-iterated. 



4. — Theory versus Experiment 



In this Section the predictions of different theoretical models for the total mesonic 
and non-mesonic weak decay rates are compared with experimental data. We mainly 
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Table I. - Sensitivity of the decay rates to the A wave function for ]^ C. 





Micr. 


Dover 


H.O. 


New 


BNL 


KEK 


KEK New 






W-S 




W-S 


[37] 


[38] 


[39-41] 




0.25 


0.25 


0.26 


0.25 


0.11 ±0.27 


0.36 ±0.13 


0.31 ±0.07 


Ti 


0.69 


0.77 


0.78 


0.82 








T2 


0.13 


0.15 


0.15 


0.16 










0.81 


0.92 


0.93 


0.98 


1.14 ±0.20 


0.89 ±0.18 


0.83 ±0.09 


Tt 


1.06 


1.17 


1.19 


1.23 


1.25 ±0.18 


1.25 ±0.18 


1.14 ±0.08 



refer to recent works. For a more specialized discussion of theoretical and experimental 
results see Ref. [5]. 

4'1. Phenom,enological a,pproach . Wc illustrate here the results which have been 
obtained for the decay rates by employing the PPM in the phenomenological approach. 

In order to evaluate the widths from Eqs. (23), (24), one needs to specify the wave 
function for the A. In Ref. [27] it has been obtained from a Woods-Saxon A-nucleus 
potential with fixed diffuseness a = 0.6 fm and with radius and depth such that it exactly 
reproduces the first two single particle eigenvalues (s and p A-levels) of the hypernucleus 
under analysis. 

A crucial ingredient of the calculation is the short range part of the strong NN 
and AN interactions, Vl, Vt, S, Pl and Pt entering Eq. (12). They are expressed by 
the functions gL,Tiq) and g^Til) reported in the appendix of Ref. [5] and contain the 
Landau parameters g' and 5^, respectively. No experimental information is available 
on g^, while many constraints have been set on g', for example by the well known 
quenching of the Gamow-Teller resonance. Realistic values of g' within the framework 
of the ring approximation are in the range 0.6-j-0.7 [22]. However, in the present context g' 
correlates not only p-h pairs but also p-h with 2p-2h states. Accordingly, the correlation 
parameters have been fixed to g' = 0.8 and = 0.4 [27] in order to reproduce the non- 
mesonic width measured for j^^C. 

Using these values for the Landau parameters, we illustrate in Table I the sensitivity of 
the calculation to the A wave hmction in j^C. In addition to the Woods Saxon potential 
(New W-S) that reproduces the s and p A-levels, other choices have been introduced: an 
harmonic oscillator wave function (H.O.) with an "empirical" frequency u), obtained from 
the s — p energy shift, the Woods-Saxon wave function of Ref. [35] (Dovc;r W S) and the 
microscopic wave function (Micr.) calculated, in Ref. [36], from a non-local self-energy 
using a realistic AN interaction. The results (in units of the free A width) are compared 
with the experimental data from BNL and KEK. By construction, the chosen g' and 
g'j^ reproduce the experimental non-mesonic width using the W-S wave function which 
gives the right s and p hyperon levels in \^C. We note that it is possible to generate 
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Table II. Mass dependence of the hypernuclear weak decay rotes (taken from Ref. [5]). 



A+ly 

A ^ 






r2 




iHe 

A ^ 

|Si 

A 

fFe 

89v 
A 


0.60 

0.25 

0.07 

0.03 

0.01 
6 X 10-3 
6 X 10-3 
1 X 10-^ 


0.27 
0.82 
1.02 
1.05 
1.12 
1.16 
1.14 
1.21 


0.04 
0.16 
0.21 
0.21 

0.21 
0.22 
0.18 
0.19 


0.91 

1.23 
1.30 
1.29 
1.35 
1.38 
1.33 
1.40 



the microscopic wave function of Ref. [36] for carbon via a local hyperon-nucleus W-S 
potential with radius R = 2.92 fm and depth Vq = —23 MeV. Although this potential 
reproduces fairly well the experimental s-level for the A in ]^G, it does not reproduce 
the p-level. A completely phenomenological A-nucleus potential, that can easily be 
extended to heavier nuclei and reproduces the experimental A single particle levels as 
well as possible, has been preferably adopted in Ref. [27]: the potential parameters 
obtained for carbon are R = 2.27 fm and Vq = —32 MeV. 

To analyze the results of Table I, we note that the microscopic wave function is sub- 
stantially more extended than all the other wave functions. The Dover's parameters [35], 
namely R = 2.71 fm and Vq = —28 MeV, give rise to a A wave function that is somewhat 
more extended than the New W-S one but is very similar to the one obtained from a 
harmonic oscillator with an empirical frequency tko = 10.9 MeV. Consequently, the non- 
mesonic width from the Dover's wave function is very similar to the one obtained from 
the harmonic oscillator and slightly smaller than the new W-S one. The microscopic 
wave-function predicts the smallest non-mesonic widths due to the more extended A 
wave-function, which explores regions of lower density, where the probability of inter- 
acting with one or more nucleons is smaller. From Table I one can also see that, against 
intuition, the mesonic width is quite insensitive to the A wave function: the different 
choices give rise to total decay widths which may differ at most by 15%. 

Using the New W-S wave fimctions and the Landau parameters g' = 0.8 and g'jy = 0.4, 
in Refs. [5,27] the calculation has been extended to hypernuclei from iHe to i°®Pb. 
Notice that, in order to reproduce the experimental s and p levels for the hyperon in 
the different nuclei one must use potentials with nearly constant depth, around 28 -r- 32 
MeV, in all but the lightest hypernucleus (iHe). For this hypernucleus the A-nucleus 
mean potential has a repulsive core and the most convenient A wave function turns out 
to be the one derived in Ref. [42] within a quark model description of iHe. 

Table II shows the resulting hypernuclear decay rates in units of the free A width. We 
observe that the mesonic rate rapidly vanishes by increasing the nuclear mass number A. 
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This is well known and it is related to the decreasing phase space allowed for the mesonic 
channel, and to smaller overlaps between the A wave function and the nuclear surface, 
as A increases. In Fig. 8 the results of Ref. [5, 27] for Fm are compared with the ones of 
Nieves-Osct [8] and Motoba-Itonaga-Bando [9,43], which were obtained within a shell 
model framework. Also the central values of the available experimental data are shown. 
Although the WFM is more reliable than the LDA for the evaluation of the mesonic 
rates (since the small energies involved in the decay amplify the effects of the nuclear 
shell structure), one sees that the LDA calculation of Ref. [27] fairly agrees with the WFM 
ones and with the data. In particular. Table II shows that the results for and ^^Si are 
in agreement with the very recent KEK measurements [41]: rM(A^C)/F^®^ = 0.31 ±0.07, 
F„- (a Si)/F^^°° = 0.046±0.011. The rcsuhs for f Ca, f Fc and fY are in agreement with 
the old emulsion data (quoted in Ref. [12]), which indicates r^-/rNM — (0.5-^ 1) x 10~^ 
in the region 40 < A < 100. Moreover, the recent KEK experiments [41] obtained the 
limit: F^- (^6Fc)/F^^°° < 0.015 (90% CL). It is worth noticing, in figure 8, the rather 
pronounced oscillations of Fm in the calculation of Refs. [9, 43], which are caused by shell 
effects. 

As a final comment to Table II, we note that, with the exception of ^He, the two- 
body induced decay is rather independent of the hypernuclear dimension and it is about 
15% of the total width. Previous works [26,44] gave more emphasis to this new channel, 
without, however, reproducing the experimental non-mesonic rates. The total width 
does not change much with A, as it is also shown by the experiment. 

In Fig. 9 the results of Table II are compared with recent (after 1990) experimental 
data for Fnm and Ft. The theoretical results are in good agreement with the data 
over the whole hypernuclear mass range explored. The saturation of the AN — > nN 
interaction in nuclei is well reproduced. 

4'2. Microscopic approach. - The results discussed in this subsection have been ob- 
tained by applying the one boson-loop formalism developed in Section 3'3.2 for nuclear 
matter. Although, in principle, one could extend this calculation to finite nuclei through 
the local density approximation, in practice this would require prohibitive computing 
times. Indeed, the latter are already quite conspicuous for the evaluation of the dia- 
grams of Fig. 7 at fixed Fermi momentum. Hence, in order to compare the results with 
the experimental data in finite nuclei, different Fermi momenta have been employed in 
the nuclear matter calculation. The average, fixed Fermi momentum, which is appropri- 
ate for the present purposes, has been obtained by weighting each local Fermi momentum 
kp (r) with the probability density of the hyperon in the considered nucleus: 



It is then possible to classify the hypernuclei for which experimental data on the non- 
mesonic decay rate are available into three mass regions (medium-light: A ~ 10; medium: 

A ~ 30-^ 60; and heavy hypernuclei: A > 200), as shown in Table III. The experimental 
bands include values of the non-mesonic widths which are compatible with the quoted 
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Fig. 8. - Mesonic width as a function of the nuclear mass number A. The results of Rof [5, 
27] (thick solid line) are compared with the calculations of Nieves-Oset [8] (dashed line) and 
Motoba-Itonaga-Bando [9,43] (solid line). Available experimental data [37,38,41] are also 
shown. See text for details on data (taJsen from Ref. [5]). 



experiments. For medium and heavy hypernuclei the available experimental data actually 
refer to the total decay rate. However, the mesonic width for medium hypernuclei is at 
most 5% of the total width and rapidly decreases as A increases. Hence, one can safely 
approximate with Fl^f^ for medium and heavy systems. In the third column of 

Table III we report the average Fermi momenta obtained with Eq. (44). According 
to the Table, in Ref. [28] the following average Fermi momenta have been employed: 
kp = 1-1 fm~^ for medium-light, kp = 1-2 fm~^ for medium and kp = 1-36 fm~^ for 
heavy hypernuclei. 

In addition to {kp), other parameters enter into the microscopic calculation of hy- 
pernuclear decay widths, which are specifically related to the baryon-meson vertices and 

to the short range correlations. With the exception of the Landau parameters g' and 
g'^, the values of these parameters have been kept fixed on the basis of the existing 
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Table III. - Average Fermi momenta for three representative mass regions. The experi- 
mental data are in units of the free A decay rate (taken from Ref. [28]). 



Medium Light: ^^B - 


0.94^ 


1.07 [37,38] 


1.08 


Medium : f Si - ^^Fe 


1.20 


^ 1.30 [39] 


~ 1.2 


Heavy: io^Bi - f^\J 


1.45-=- 


1.70 [45,46] 


1.36 
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1.1 1.2 1.3 1.4 

kjr (fm~^) 

Fig. 10. - Dependence of the non-mesonic width on the Fermi momentum of nuclear matter. 

The sohd curves refer to the one-boson-loop approximation (with g' — 0.7, 0.8, 0.9 from the top 
to the bottom), while the dashed lines refer to the ring approximation (g' = 0.5,0.6,0.7). The 
experimental data are also shown (taken from Ref. [28]). 



phenomenology (for example in the analysis of quasi-elastic electron-nucleus scattering, 
spin-isospin nuclear response functions, etc). For the complete list of these quantities we 
refer the reader to Ref. [28]. 

Instead, the zero energy and momentum limits of the strong NN and AA'' correlations, 
g' and g^, are considered as free parameters. Wc remind the reader that the physical 
meaning of these Landau parameters is different in the present scheme (see Fig. 7) with 
respect to the customary phenomenology based on the ring approximation. 

Fixing g'j^ = 0.4, in ring approximation the experimental decay rates can be repro- 
duced by using g' values which are compatible with the existing literature (0.5 < g' < 
0.7). However, larger g' values (0.7 < g' < 0.9) appear to be more appropriate in the 
framework of the full OBL approximation. In addition, the OBL calculation allows for 
a good description (keeping the same g' value) of the rates in the whole range of kp 
considered here. 

In Fig. 10 we see the dependence of the non-mesonic widths on the Fermi momentum. 
The solid lines correspond to the OBL approximation, with g' = 0.7,0.8,0.9 from the 
top to the bottom, while the dashed lines refer to the ring approximation, with g' = 
0.5,0.6,0.7, again from the top to the bottom. One can then conclude that for the 
OBL calculation the best choice for the Landau parameters is the following one: g' = 
0.8, g'jy = 0.4. Accidentally, this parameterization turns out to be the same employed 
in the phenomenological evaluation of the 2p-2h polarization propagator, although the 
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Table IV. 



Mesonic decay rate for \He. 



Ref. 



Tee 
A 



Model 



Oset-Salcedo 1985 [24] 0.65 

Oset-Salccdo-Usmani 1986 [49] 0.54 

Itonaga-Motoba-Bando 1988 [43] 0.331 ^ 0.472 

Motoba et al. 1991 [47] 0.608 

Motoba 1992 [48] 0.61 

Straub et al. 1993 [42] 0.670 

Kumagai-Fuse et al. 1995 [50] 0.60 



PPM 

PPM 

WFM 
WFM + Quark Model 

WFM 
WFM + Quark Model 

WFM 



Exp BNL 1991 [37] 
Exp KEK 2004 [51] 



0.59t°it 
0.541 ± 0.019 



role of strong short range correlations is different in the two approaches. 

4'3. Other evaluations. - In this subsection we briefly discuss how other calculations 
— ^performed within the Polarization Propagator Model (PPM) and the Wave Function 
Model (WFM) previously illustrated — compare with existing data. 

In Table IV results for the mesonic decay width for ^He are reported. All the theoreti- 
cal evaluations agree with the BNL datum but discrepancies appear if one compares them 
with the more accurate KEK observation. A repulsive core in the A-a mean potential 
(used in all but the calculation of Ref. [24]) is favoured. Moreover, it comes out naturally 
in the quark model descriptions of Refs. [42,47]. The results of Refs. [43,48] refer to the 
use of different pion-nucleus optical potentials. On this point, we note that the more 
recent datum is able to discriminate between different pion-nucleus optical potentials. 
Thanks to such precise determinations of Fm, novel information on this potential will 
be possibly extracted in the next future. The experimental and theoretical values of 
the ratio r^-/r^o, not shown in the table, do not deviate much from the Al =1/2 
value (= 2) for free decays. We expect this result, since ^He has a closed shell core with 
N = Z. 

The theoretical results for the mesonic rate for j^^C, reported in Table V, are fairly 

compatible with data. The most notable exception is the calculation of Ref. [52], supply- 
ing a decay rate which strongly underestimates the KEK data. Also in this case, accurate 
measurements such as the one of Ref. [51] would provide important information on the 
pion-nucleon optical potential. The estimates obtained with the WFM of Refs. [8,9, 43] 
are consistent with the experimental ratio (F^o /T^- )^^p ~ 1 ^ 2 > (F^o /F^- Y''"'' = 1/2, 
which reflects the particular nuclear shell structure of j^C. 

Table VI refers to the non-mesonic rate for ^He, j^^C and nuclear matter. The dif- 
ferent AA'' nN potentials used in the calculations reflect in a broad spectrum of 
predictions. In particular, Dalitz [55] (and more recently Jun et al. [56]) used a phe- 
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Table V. - Mesonic decay rate for ]^ C. 



Ref. 




Model 


Oset-Salcedo 1985 [24] 


0.41 


PPM 


Itonaga Motoba-Bando 1988 [43] 


0.233^0.303 


WFM 


Ericson-Bando 1990 [53] 


0.229 


WFM 


Nieves-Oset 1993 [8] 


0.245 


WFM 


Itonaga-Motoba 1994 [9] 


0.228 


WFM 


Ramos et al. 1995 [26] 


0.31 


PPM 


Zhou-Piekarewicz 1999 [52] 


0.112 


Relativistic PPM 


Albertus et al. 2003 [54] 


0.25 


WFM 


Exp BNL 1991 [37] 


0.11 ±0.27 




Exp KEK 1995 [38] 


0.36 ±0.13 




Exp KEK 2004 [41] 


0.313 ±0.070 




Exp KEK 2004 [51] 


0.288 ±0.017 





nomenological model in which the OPE model at large distances is supplemented by a 
4 baryon point interaction (4BPI) for the short range AiV nN interactions. Ref. [20] 
used a hybrid model in which the long range interactions are treated in terms of the 
OPE, while the short range interactions are described by a 6-quark cluster model which 
includes both A/ = 1/2 and Al = 3/2 components. The PPM of Ref. [24] overesti- 
mate the data, although the decay rate is reduced when a more realistic A-wave function 
(less superimposed with the nuclear core) is used [27, 28, 49], and particularly if different 
short range correlations are considered [27,28] (see results discussed in subsections 4'1 
and 4'2). Antisymmetrization of the final nucleons, as in Ref. [18], would also moderately 
decrease the non-mesonic rate. Sasaki et al. [21] treated the non-mesonic decay within 
a direct quark (DQ) model combined with a OME potential containing pion and kaon 
exchange. In their model the NN and AA^ repulsion at short distance originates from 
quark exchange between baryons (induced by the quark anti-symmetrization) and gluon 
exchange between quarks. The intervals shown for the OME calculation of Ref. [17] 
correspond to the use of different Nijmegen models for the hadronic coupling constants. 
Itonaga et al. [16] considered the correlated two-pion-exchange. 

Large part of the predictions of Table VI agree with the available data. The non- 
mesonic width in j^^C seems to be reduced by a factor of about 2 with respect to the 
nuclear matter value. Prom inspection of data on heavy hypernuclei one concludes that 

realistic values of the A decay rate in nuclear matter lie in the range 1.5 -j- 2. From both 
experiment and theory, it also follow that rNM(A^C) ~ 2rNM(AHe). 
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Table VI. - Non-mesonic decay rate for %He, C and nuclear matter in units of the free 
A decay width. 



Ref. and Model 


iHe 


A ^ 


Nuclear Matter 


T^^l'j- 1 ATO feci 

Dalitz 1973 [55J 


0.5 




2 


^WriVl: yjrtj + 4r5r^lJ 








Cheung et al. 1983 [20J 




1.28 


3 


(WFM: hybrid) 








Oset-Salcedo 1985 [24] 


1.15 


1.5 


2.2 


(PPM: Correlated OPE) 








Oset-Salcedo-Usmani 1986 [49] 


0.54 






(PPM: Correlated OPE) 








Sasaki et al. 2000 [21] 


0.519 




2.456 


(WFM: TT + DQ) 








Jun et al. 2001 [56] 


0.426 


1.174 




(WFM: OPE + 4BPI) 








Jido et al. 2001 [18] 




0.769 




(PPM: TT + K + 27r + w) 








Parreno-Ramos 2001 [17] 


0.317^0.425 


0.554 ^ 0.726 




(WFM: ■K + p + K + K*+uj + 7]) 








Itonaga et al. 2002 [16] 


0.422 


1.060 




(WFM: TT + 2tt/p + 27r/cr + uj) 








Exp BNL 1991 [37] 


0.41 ±0.14 


1.14 ±0.20 




Exp CERN 1993 [45] 






p+Bi: 1.46tjj3 
p+U: 2.021 J:^^ 


Exp KEK 1995 [38] 




0.89 ±0.18 




Exp KEK 1995 [57] 


0.50 ±0.07 






Exp COSY 1998 [46] 






p+Bi: l.&^taZ 


Exp KEK 2000 [39,40] 




0.83 ±0.11 


fYc: 1.22 ±0.08 


Exp COSY 2001 [58] 






p±Au: 2.02tEJ-^6 


Exp COSY 2001 [59] 






p+\5: i.9i;;;-2« 


Exp KEK 2004 [-ll] 


0.406 ± 0.020 


0.9.^)3 ± 0.0:-!2 





5. - The ratio r„/rp 

Up to very recent times, the main challenge of hypernuclear weak decay studies has 
been to provide a theoretical explanation of the large experimental values for the ratio 

Fn/Fp between the neutron and the proton-induced decay widths. Until recently, the 
large uncertainties involved in the extraction of the ratio from data did not allow to 
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reach any definitive conclusion. The data were quite hmited and not precise due to the 
difficulty of detecting the products of the non-mesonic decays, especially the neutrons. 
Moreover, up to now it has not been possible to distinguish between nucleons produced 
by the one-body induced and the (non- negligible) two body induced decay mechanism. 
However, due to recent theoretical [16-18,21,60] and experimental [61-63] progress, we 
are now towards a solution of the r^/Fp puzzle. In this section we review these important 
developments. 

The one-pion-exchange approximation with the AI =1/2 rule provides small ratios: 



(45) 



p nOPE 



0.05 -=-0.20, 



for all the considered systems. This is due to the A7 =1/2 rule, which fixes the vertex 

ratio Vf^^-p/V\Tr"n = ~\/2 (both in S- and P-wave interactions), and to the particular 
form of the OPE potential, which has a strong tensor component requiring isospin np 
pairs in the antisymmetric final state. In p-shell and heavier hypernuclei the relative 
AN L = 1 state is found to give only a small contribution to tensor transitions for the 
neutron-induced decay, so it cannot improve the OPE ratio. The contribution of the AA^ 
L = 1 relative state to Fnm seems to be of about 5 15% in p-shell hypernuclei [17]. 
For these systems we expect the dominance of the S-wave interaction in the initial state, 
due to the small AA'^ relative momentum. By using a simple argument about the isospin 
structure of the transition AA^ nN in OPE, it is possible to estimate that for pure 
A7 = 3/2 transitions (V^^— p/VA,rOra = 1/v^) the OPE ratio can increase up to about 
0.5. However, the OPE model with AI =1/2 couplings has been able to reproduce the 
one-body stimulated non-mesonic rates Fi = F„ -f Tp for light and medium hypernuclei 
[16-18,21,27,28]. 

Other interaction mechanisms beyond the OPE might then be responsible for the 

overestimation of Tp and the underestimation of r„. Many attempts have been made up 
to now in order to solve the F„/Fp puzzle. We recall here the inclusion in the AA^ nN 
transition potential of mesons heavier than the pion [14-18] , the inclusion of interaction 
terms that explicitly violate the AI =1/2 rule [19] and the description of the short range 
baryon-baryon interaction in terms of quark degrees of freedom [20, 21], which automat- 
ically introduces AI = 3/2 contributions. It is important to note that a few calculations 
with A A'' — > nN transition potentials including heavy-meson-exchange and/or direct 
quark contributions [16-18,21] have improved the situation, without providing, never- 
theless, a satisfactory explanation of the origin of the puzzle. The tensor component 
of ii'-exchange has opposite sign with respect to the one for 7r-exchange, resulting in a 
reduction of Tp. The parity violating AN{^Si) — > nN{^Pi) transition, which contributes 
to both the n- and p-induced processes, is considerably enhanced by JT-exchange and 
direct quark mechanism and tends to increase F„/Fp. 

In table VII we summarize the calculations that predicted ratios considerably en- 
hanced with respect to the OPE values. Experimental data are given for comparison. 
Almost all calculations reproduce the observed non-mesonic widths F„ -|- Fp, as one can 
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A ^ 



0.701 

0.343 ^ 0.457 
0.386 



0.53 
0.288^0.341 
0.368 



0.93 ±0.55 
1.97 ±0.67 
0.40 ±0.11 



1 00 + 1.12 
1.00_o.81 

1.871?:?^ 



0.87 ±0.23 
0.38 ±0.14 



Rcf. and Model 



Sasaki et al. 2000 [21] 
TT + K+BQ 
Jido et al. 2001 [18] 

TT + K + 2TT + UJ 

Parreiio-Ramos 2001 [17] 
TT + p + K + K*+u; + r] 
Itonaga et al. 2002 [16] 
TT ± 27r//9 ± 27r/cr ± w 



BNL 1991 [37] 
KEK 1995 [38] 
KEK 1995 [57] 
KEK 2004 [41] 
KEK 2004 (coincidence) [61] 
(with Garbarino et al. [60, 64] analysis) 



see in Table VIII [we remind the reader that the experimental data should also include 
(at least a part) of the two-nucleon induced decay rate]. Although no calculation is able 
to explain the old data on Fn/Fp, some predictions are in agreement with the recent 
determinations [60,64] from KEK nucleon coincidence data [61]. We discuss in detail 
these analyses in 5 '2. 2. 

5'1. Two-body induced decay and nucleon final state interactions. - The analysis of 
the ratio F^/Fp is influenced by the two-nucleon induced process ANN — > nNN, whose 
experimental identification is rather difficult and it is a challenge for the future. By 
assuming that the meson produced in the weak vertex is mainly absorbed by an isoscalar 
np correlated pair (quasi-deuteron approximation), the three-body process turns out to 
be Anp — > nnp, so that a considerable fraction of the measured neutrons could come 
from this channel and not only from An — > nn and Ap — > np. In this way it might 
be possible to explain the large experimental F„/Fp ratios, which originally have been 
analyzed without taking into account the two-body stimulated process. Nevertheless, the 
situation is far from being clear and simple, both from the theoretical and experimental 
viewpoints. The new non-mesonic mode was introduced in Ref. [44] and its calculation 
was improved in Ref. [26] , where the authors found that the inclusion of the new channel 
would bring to extract from the experiment even larger values for the F„/Fp ratios, thus 
worsening the disagreement with the theoretical estimates. However, in the hypothesis 
that only two out of the three nucleons coming from Anp nnp are detected [65] , the 
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Table VIII. - Non-mesonic width r„ + Tp (in units o/r^**®^. 







A ^ 


Sasaki et al. 2000 [21] 


0.519 










Jido et al. 2001 [18] 




0.769 


71 + K + 2TT + UJ 






Parreno-Ramos 2001 [17] 


0.317^0.425 


0.554^0.726 


'K + p + K + K*+uj + 'n 






Itonaga et al. 2002 [16] 


0.422 


1.060 


7r + 27r/p + 27r/o- + w 






BNL 1991 [37] 


0.41 ±0.14 


1.14 ±0.20 


KEK 1995 [38] 




0.89 ±0.18 


KEK 1995 [57] 


0.50 ±0.07 




KEK 2000 [40] 




0.83 ±0.11 


KEK 2004 [51] 


0.406 ± 0.020 


0.953 ± 0.032 



reanalysis of the experimental data would lead back to smaller ratios. 

These observations show that r„/rp is sensitive to the detailed kinematics of the non- 
mesonic processes and to the experimental threshold for nucleon detection. In Ref. [66] 
the nucleon energy distributions have been calculated by using a Monte Carlo simulation 
to describe the nucleon rescattcring inside the nucleus: the ratio Fn/Fp has been taken as 
a free parameter and extracted by comparing the simulated spectra with data. The mo- 
mentum distributions of the primary nucleons were determined within the polarization 
propagator scheme. In their way out of the nucleus, the nucleons, due to the collisions 
with other nucleons, continuously change energy, direction, charge, and secondary nucle- 
ons are emitted as well. Then, the energy distribution of the observable nucleons, which 
also loose their energy by the interactions with the experimental set-up, is quite different 
from the one at the level of the primary nucleons. The nucleons from two-body stimu- 
lated decays appear mainly at low energies, while those from the one-nucleon stimulated 
process peak around 75 MeV only for light hypernuclei such as ^He. 

Quite recently, at KEK-E307 [67], the proton spectra for j^^C, f^Si and ^^Fe were 
measured and compared with theoretical simulations of the intranuclear cascades after 
the weak processes, obtained with the MC code of Ref. [66] . Corrections for the detector 
geometry and the nucleon interactions inside the target and detector materials were 

necessary: they were implemented through a GEANT MC code. After fitting the KEK- 
E307 spectra using F„/Fp as a free parameter, the following results were obtained [41] 
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by neglecting the two-nucleon induced decay channel: 



(46) 




0.87 ±0.23 





All the existing calculations underestimate the ^^C result (see Table VII). The Tn/Tp 
ratios of Eq. (46) confirm the results of previous experiments: the neutron- and proton- 
induced decay rates appear to be of the same order of magnitude over a large hypernuclear 
mass number range. However, since the new experiments have significantly improved the 
quality of the data, small values of Tn/Tp (say smaller than 0.5 for j^^C, as predicted by 
theory) seem to be excluded. 

5'2. Towards a solution of the Tn/Tp puzzle. - Fortunately, recent important devel- 
opments have contributed to approach the solution of the Fn/Fp puzzle. This progress 
has been based on the following main ideas: 

1) The proton spectra originating from neutron- and proton-induced processes (and, 
eventually, from two-nucleon stimulated decays) are added incoherently in the 
Monte Carlo intranuclear cascade code used to determine F„/Fp from data. In 
this way a possible quantum-mechanical interference effect between the two chan- 
nels is lost. Therefore, extracting F„/Fp from data with the help of a classical 
intranuclear cascade calculation may not be a clean task. The consequences of this 
idea have been explored in Ref. [60]. 

2) In order to perform, as desirable, a direct measurement of F„/Fp one needs to 
detect the outgoing neutrons. In principle, neutron spectra can be measured down 
to about 10 MeV kinetic energy since they are less affected than the proton ones 
by energy losses in the target and detector materials. Besides, the joint observa- 
tion of proton and neutron spectra could help to disentangle the set-up material 
effects from the nucleon FSI occurring inside the residual nucleus. A recent ex- 
periment, KEK E369, measured neutron spectra from and ^^Y non-mesonic 
decays [62]. With other experiments, KEK E462 and KEK E408 [63], proton and 
neutron spectra could be simultaneously measured for ^He and \^C. An analysis 
of these data has been done in Ref. [60] and is discussed in 5 '2.1. 

3) One could avoid the possible deficiencies of the single nucleon spectra measure- 
ments discussed in point 1) by employing nucleon-nucleon coincidence measure- 
ments, which should indeed permit a cleaner extraction of F„/Fp. Coincidence 
observations are also expected to be less affected by FSI effects and to help in the 
direct observation of two-nucleon induced decay events. In the experiments KEK- 
E462 and KEK-E508 [61], nn and np angular and energy correlation measurements 



40 



W. M. Alberico and G. Garbarino 



have been performed for the decay of ^Hc and ]^C. Other experiments, at BNL 
[68] and J-PARC [69], will determine Tn/Tp for s-shell hypernuclei, again by nn 
and np coincidence measurements. More data are expected in the near future from 
Da$ne [70]. 

5 '2.1. Analyses of single nucleon spectra. In Ref. [60] a OME model for the AA^ nN 
transition in finite nuclei was incorporated in the intranuclear cascade code of Ref. [66] 
for the calculation of single and coincidence nucleon distributions from hypernuclear non- 
mesonic weak decay. The OME weak transition potential [17] contained the exchange 
of TT, /?, K, K* , Lu and r) mesons. Strong couplings and strong FSI acting between 
the weak decay nucleons were taken into account by means of a nN wave function 
from the Lippmann-Schwinger equation obtained with NSC97 (versions "a" and "f) 
potentials [71]. The two-nucleon stimulated channel is also taken into account, using the 
polarization propagator method and treating the nuclear finite size effects by means of 
a local density approximation. 

Single neutron and proton spectra have been calculated and compared with data for 
^He and ]^C. In Fig. 11 we show results of Ref. [60] for the kinetic energy neutron dis- 
tribution from based on two models [OPE and OMEf (using the NSC97f potential)], 
which predict quite different Tn/Tp ratios. The KEK-E369 spectrum [62] is well repro- 
duced by the calculations. Unfortunately, the dependence of the neutron spectrum A„ 
on variations of Tn/Tp is very weak (the same is true also for the proton spectrum) and 
a precise extraction of the ratio from the KEK-E369 distribution is not possible. The 
little sensitivity of A„ and Np to Tn/Tp is mainly due to the fact that these numbers 
are normalized to the same total non-mesonic decay rate (i.e., per non-mesonic weak 
decay). The non-normalized nucleon spectra, Sn = A^Fnm and Sp = ApFNM [see 
Eq. (47)], have indeed a stronger dependence on Tn/Tp. As a consequence, in order to 
discriminate between different weak decay models, one should separately compare the 
complementary observable, Fnmj with experiment. For ^^C, the calculations of Ref. [60] 
supply Fnm = F„ + Fp + F2 = 0.91 or 0.69 when model OMEa or OMEf is used, which 
agree quite well with the experimental determinations (see Table VIII). 

Let us now introduce the number of nucleons of the type N (N = n,p) produced 
in n-induced (A^^"), p-induced {Nj^^} and two-nucleon induced (A|P) decays. If we 
normalize these quantities per n-induced, p-induced, and two-nucleon induced decay, 
respectively, the total number of nucleons of the type N (normalized per non-mesonic 
weak decay) is given by: 



(47) Nn = 



Aip" r„ + NjfP Tp + NfP T2 



Tn ~t- Fp 

AjKn^nn , j^^rAp^np , AjAnp^nnp 



where A^"^"", A^^^"^ and 7v^"p^""p have obvious meaning and are shown in Fig. 11 

for the neutron spectrum from }^C. All these nucleon numbers can be considered either 
as being functions of the nucleon kinetic energy, Nn{Tn), as it is in Fig. 11, or as 
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total OPE 

from n-induced 

from p-induced 

from 2N-induced 

• KEK-E369 




T„ (MeV) 

Fig. 11. - Single neutron kinetic energy spectra for the non-mesonic weak decay of ]^C. The total 
spectrum iV„ (normalized per non-mesonic weak decay) has been decomposed in its components 
iV4"^"", N^^^"" and iV^"''^""^ according to Eq. (47) (taken from Ref [60]). 



0.5 
0.4 



the corresponding integrated quantities, Nn = J dTNNN{TN)- By construction, Nj:^'^, 
A^^^P and iV2B (ivAn-"„^ N^p^'^p and ArA"P-»"P) depend {do depend) on the 

interaction model employed to describe the weak decay. 

The problem of the small sensitivity of iV„ and Np to variations of Tn/Tp can be 
overcome if one concentrates on another single nucleon observable, the ratio Nn/Np. 
The number of weak decay neutrons and protons produced per non-mesonic weak decay 
are given by N:^'^ = (2r„ + rp)/rNM and N^'^ = Tp/r^M, respectively. Thus, the ratio 
Tn/Tp can be obtained as: 

Due to two-body induced decays and (especially) nucleon FSI, one expects the inequality: 
(49) I^4(^-l).«,[AT„.ATd 

to be valid in a situation, such as the experimental one, in which particular intervals of 
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Table IX. - Predictions of Ref. [60] for the quantity Ri of Eq. (4-9) for \He corresponding 
to different nucleon kinetic energy thresholds T}^ and to the OPE, OMEa and OMEf 
models. 








K (MeV) 
30 


60 




OPE 


0.04 


0.13 


0.16 


0.09 


OMEa 


0.15 


0.32 


0.39 


0.34 


OMEf 


0.19 


0.40 


0.49 


0.46 


KEK-E462 [63] 






0.59 ±0.11 





variability of the neutron and proton kinetic energy, AT„ and ATp, are employed in the 
determination of the observable numbers iV„ and Np. 

The results of Ref. [60] clearly show a pronounced dependence of i?i on AT„ and 
ATp (sec Tabic IX); Nn/Np turns out to be much less sensitive to FSI effects and varia- 
tions of the energy cuts than Nn and Np separately. The OMEf prediction of Table IX 
corresponding to = 60 MeV is compatible with the very recent determination by 
KEK-E462 [63]. 

5'2.2. Analyses of coincidence spectra. Due to the reduction of interferences and FSI 
effects, nucleon correlation analyses are expected to provide a cleaner determination of 
r„/rp than single nucleon observables [5]. Ref. [60] evaluated double-nucleon energy 
and angular correlations and analyzed the data obtained by the experiments KEK-E462 
and KEK-E508 [61] for ^He and j^C. 

In Fig. 12 we report the prediction of Ref. [60] for the kinetic energy correlation of np 
pairs emitted in the non-mesonic decay of ]^C. To facilitate a comparison with experi- 
ments, whose kinetic energy threshold for proton (neutron) detection is typically of about 
30 MeV (10 MeV), and to avoid a possible non-realistic behaviour of the intranuclear 
cascade simulation for low nucleon energies, the distributions correspond to T„, Tp > 30 
MeV. A narrow peak, mainly originated by the back-to-back kinematics (cos 9np < —0.8) 
of the one-nucleon induced decay, is predicted around 155 McV, i.e., close to the Q-value 
expected for the proton-induced three-body process j^C ^^"B + n + p. A broad peak, 
predominantly due to Ap np or An nn weak transitions followed by the emission of 
secondary (less energetic) nucleons, has been found around 110 MeV for cos^„p > —0.8. 

Figure 13 shows the opening angle correlations of nn, np and pp pairs from the non- 
mesonic weak decay of and for a 30 MeV energy threshold. The peaking structure at 
cos ^ ~ — 1 for the nn (np) distribution is a clear signal of the back-to-back kinematics 

of the neutron (proton) induced decay. For a discussion of the different effects of FSI in 
the energy and angular coincidence distributions for ^He and j^C see Ref. [60]. 
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Fig. 12. - Kinetic energy correlations of np pairs emitted per non-mesonic decay of (taken 
from Ref. [60]). 



The ratio F^/Fp is defined as the ratio between the number of weak decay nn and np 
pairs, A^^^ and N^^- However, due to two-body induced decays and (especially) nucleon 
FSI effects, one has: 

r N""'^ N 

(50) ^^j^^^^R2 [M,,, Ar„, AT,] 

when the observable numbers Nnn and Nnp are determined by employing particular 
intervals of variability of the pair opening angle, A9i2, and the nucleon kinetic energies, 
AT„ and ATp. The discussion of Rcf. [60] proves that Nnn /Nnp is much less sensitive to 
FSI effects and variations of the energy cuts and angular restrictions than Nnn and Nnp 
separately. 

The numbers of nucleon pairs N^n discussed up to now and normalized per non- 
mesonic weak decay are related to the corresponding quantities for the one-nucleon 
(-^Aw) two-nucleon (NfPj^) induced processes [the former (latter) being normal- 
ized per one-body (two-body) stimulated non-mesonic weak decay] via the following 
equation: 

atIB p _i_ /\/'2B 

[bL) I^NN — f7+T2 ~ ' 
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Fig. 13. - Opening angle correlations of nn, np and pp pairs emitted per non-mesonic decay of 
i^C (taken from Ref. [60]). 



where: 



(52) 



^^NN — 



lylBn p 



^IBpp 



and the remaining A'''s have obvious meaning. Therefore, the quantities NjPj^, and 



Ar2B 



(^wjv^""' ^nn*""^ ^nd N^^^™^) do not depend {do depend) on the interaction 



NN ' ^"NN '^^^^ ^"NN 

model employed to describe the weak decay. 

In Table X the ratio 7V„„/iV„p predicted by the OPE, OMEa and OMEf models of 
Refs. [60, 64] for ^He and j^C is given for a nucleon energy threshold of 30 MeV and for 
the back-to-back kinematics (cos O^n < —0.8). The predictions of the different weak 
decay models for F^/Fp are also quoted. The results of the OMEa and OMEf models are 
in agreement with the preliminary KEK-E462 and KEK-E508 data [61]: this comparison 
provides an indication for a ratio Tn/Tp ~ 0.3 in both hypernuclei. 

A weak-decay-model independent analysis of KEK coincidence data has been per- 
formed in Ref. [60]. The six weak-interaction-model independent quantities iV^)^", N^n^, 
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Table X. - Predictions of Refs. [60, 64] for the ratio Nnn/Nnp for \He and \^ C. They 
correspond to an energy thresholds of 30 MeV and to the back-to-back kinematics 
(cos Onn < -0.8). The (preliminary) data are from KEK-E462 and KEK-E508 [61]. 





iHe 


r /r 


12 p 

^nnf-^np 




OPE 


0.25 


0.09 


0.24 


0.08 


OMEa 


0.51 


0.34 


0.39 


0.29 


OMEf 


0.61 


0.46 


0.43 


0.34 


EXP 


0.45 ±0.11 




0.40 ±0.09 





iV^B, iYi^"^, N^^P and N^^ of Eqs. (51) and (52) are used to evaluate Tn/Tp as: 

^IBp , ^2b£2 _ f ^iBp ^2b£2\ N„ 



p nn ' nn p I np ' np 'p j ]\t 
^^2^ _ 1 \ ^ 1 / ^^np 

( TW-lBn , TUSB^S \ ^nn _ ArlBn _ /v/-2b£2 



1 / -''np 1 

[which can be easily obtained from Eqs. (51) and (52)], using appropriate values for 
r2/ri and data for Nnn/Nnp. For the KEK values of Table X the following ratios have 
been determined: 

(54) ^ (iHe) = 0.40 ± 0.11 if Fa = , ^ (iHe) = 0.27 ± 0.11 if ^ = 0.2, 



(55) ^ (a C) = 0.38 ± 0.14 if Fj = , ^ (i^C) = 0.29 ± 0.14 if = 0.25, 
1 p 1 p 1 1 

where the non-vanishing values adopted for r2/ri are predictions obtained in Ref. [27] 
within the polarization propagator method in local density approximation. The Tn/Tp 
values of Eqs. (54) and (55) are in agreement with the pure theoretical predictions of 
Rcfs. [16 18] but arc substantially smaller than those obtained experimentally from sin- 
gle nucleon spectra analyses (see Table VII). Actually, all the previous experimental 
analyses of single nucleon spectra, supplemented in some cases by intranuclear cascade 
calculations, derived Fn/Fp values in disagreement with all existing theoretical predic- 
tions. The fact that the calculations of Refs. [60, 64] reproduce coincidence data for values 
of Fn/Fp as small as 0.3 -j- 0.4 could signal the existence of non-negligible interference 
effects between the n- and p-induced channels in those old single nucleon data. 

In our opinion, the achievements discussed in this subsection clearly exhibit the in- 
terest of analyses of correlation observables and represent an important progress towards 
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the solution of the r^/Fp puzzle. Forthcoming coincidence data from KEK, BNL [68], 
J-PARC [69] and FINUDA [70] could be directly compared with the results discussed 
here and in Refs. [60,64]. This will permit to achieve better determinations of Tn/Tp 
and to establish the first constraints on r2/Fi. 

6. — Hypernuclei of the s— shell and A7 =1/2 rule 

The analysis of the non-mesonic decays for s-shell hypernuclei offers an important 
and complementary tool for the solution of the F„/Fp puzzle. In addition, it is very 
suitable for testing the validity of the related A/ = 1/2 rule. 

For s-shell hypernuclei the AA^ pair is necessarily in the L = relative state, thus 
the only possible h.N — > nN transitions are the following ones (we use the spectroscopic 
notation 

(56) 



'So 


ilf = 


1) 




ilf = 


1) 


's. 


ilf = 


0) 


'Pi 


{If = 


0) 


'Pi 


(If- 


1) 




{If = 


:0) 



The An nn process has final states with isospin If = 1 only, while for Ap — > np both 

If = I and If = arc allowed. 

6'1. Phenomenological model of Block and Dalitz. ~ In the following we discuss an 
analysis performed by the authors [72] in order to explore the validity of the AJ =1/2 
rule in the one-nucleon induced A-decay. This analysis is based on a phenomenological 
model duo to Block and Dalitz [73], which we briefly outline now. 

The interaction probability of a particle which crosses an infinite homogeneous system 
of thickness ds is, classically, dP = ds/X, where A = l/{crp) is the mean free path of the 
projectile, a is the relevant cross section and p is the density of the system. Then, if we 
refer to the process AN nN, the width Fnm = dP^N^nN /dt can be written as: 

Fnm = vap, 

V = ds/dt being the A velocity in the rest frame of the homogeneous system. For a finite 
nucleus of density p{r), after introducing a local Fermi sea of nucleons one can write, 
within the semiclassical approximation: 



■ NM 



(va) / dfp{f) I V'A(r) 



where ^->\{r) is the A wave hmction in the hypernucleus and () denotes an average over 
spin and isospin states. In the above equation the nuclear density is normalized to the 
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mass number A = N + Z, hence the integral gives the average nuclcon density pA at the 
position of the A hyperon. In this scheme, the non-mesonic width Fnm — r„ + Tp of the 
hypernucleus ^'^^Z turns out to be: 

P (A+ir,. _ NRnit^'Z) + ZRpii+'Z) 
J-nmU ^) — Pa, 

where i?„ (i?p) denotes the spin- averaged rate for the neutron-induced (proton-induced) 
process appropriate for the considered hypernucleus. 

Furthermore, by introducing the rates Rnj for the spin-singlet (i?nO) Rpo) and spin- 
triplet {Rnl, Rpl) elementary AA'' — > nN interactions, the non-mesonic decay widths of 
s-shell hypernuclei are [73]: 

(57) rNM(AH) = (3i?„o + Rnl + SRpO + Rpl) 

o 

TNiiidH) = (-RnO + 3i?„i + 2Rpo) 
rNM^He) = (2i?„o + RpO + 3-Rpi) 

D 

rNM(AHe) = (-RnO + 3Rnl + RpO + ^Rpl) 

o 

These relations take into account that the total hypernuclear angular momentum is for 
and ^He and 1/2 for and %ile. In terms of the rates associated to the partial-wave 
transitions (56), the -Rjvj's of Eqs. (57) read: 

-RnO = -Rn("'^S'o) + -Rn('^-Fo), 
-RpO = -Rp(^'5'o) + -Rp(^Po), 
Rnl = Rn^Pl), 

Rpl = RpCSi) + RpQPi) + Rp^Pi) + RpCDi), 

the quantum numbers of the nN final state being reported in brackets. 

If one assumes that the AA'' nN weak interaction occurs with a change A7 = 1/2 
of the isospin, the following relations (simply derived by angular momentum coupling 
coefficients) hold among the rates for transitions to // = 1 final states: 

(58) RnCSo) = 2RpCSo), RnCPo) = 2RpCPo), Rn^Pi) = 2Rp{^Pi). 



Hence: 
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Experimental data for the non-mesonic weak decay of s -shell hypernuclei. 



XR 0.22 ± 0.09 reference value 

0.17 ±0.11 KEK [11] 

0.29 ±0.14 [73] 

^He 0.04 ±0.02 0.16 ±0.02 0.20 ±0.03 0.25 ±0.13 BNL [74] 

|He 0.20 ±0.11 0.21 ±0.07 0.41 ± 0.14 0.93 ±0.55 BNL [37] 



For pure AJ = 3/2 transitions, the factors 2 in Eqs. (58) are replaced by 1/2. Then, by 
further introducing the ratio: 

_ {If = 1\\A,/,\\I, = 1/2) 



{If = l\\A,/^\\Ii = 1/2) 

between the A/ =1/2 and A/ = 3/2 AA^ nN transition amplitudes for isospin 1 final 
states (r being real, as required by time reversal invariance), for a general A/ = 1/2-3/2 
isospin mixture one gets: 

.QQ. Rra ^ 4r^ - 4r ± 1 ^ Rno ^ 4^-2 - 4r ± 1 

^ ' Rpl 2r2 + 4r + 2 + 6A2 - iipo 2r2+4r + 2' 



where: 
(61) 



(J/=0Pi/2||/. = l/2) 

(J/ = l||A3/2||/. = l/2)- 



The partial rates of Eq. (60) supply Fn/Fp for s-shell hypernuclei through Eqs. (57). 

By using Eqs. (57) and (60) together with available data it is possible to determine 
the spin-isospin structure of the AN — > nN interaction without resorting to a detailed 
knowledge of the interaction mechanism. 

6'2. Experimental data and AI = 1/2 rule. - In Ref. [72] a phenomenological analysis 
of the data summarized in Table XI is employed to study the possible violation of the 
AI = 1/2 rule in the process AA^ nN. 

Unfortunately, no data are available on the non-mesonic decay of hypertriton (^H) 
and on F„/Fp for ^H. Indeed, we shall see in the following that future measurements of 
r„/Fp for at BNL [68] and J-PARC [69] will be of great importance for testing the 
A/ = 1/2 rule. The BNL data [37,74] for ^He and ^He of Table XI together with the 
reference value for have been used in the analysis of Ref. [72]. This last number is 
the weighted average of two previous estimates [11,73], which have not been obtained 
from direct measurements but rather by using theoretical constraints. One has then 5 
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independent data which allow to fix, from Eqs. (57), the 4 rates Rn,.j and p^. Indeed, 
the average nucleon density at the A position for ^He, also entering Eqs. (57), has been 
estimated to be p4 = 0.045 fm^'' by employing the A wave function of Ref. [42] (which 
was obtained through a quark model description of the AA^ interaction) and the gaussian 
density for ^He that reproduces the experimental mean square radius of the nucleus. For 
and XRe, instead, no realistic hyperon wave function is available and we can obtain 
the value ps = 0.026 fm~^ from the data of Table XI, by imposing that [see Eqs. (57)]: 

rp(iHc) ^ 3p4 
TpdHe) 4 PS- 

The best choice to determine the rates Rn,j by fitting measured values corresponds 
to use the relations for the observables: 

rNM(AH), rNM(AHe), rNM(AHe), — (^He), 

which have the smallest experimental uncertainties. After solving these equations we 

obtained the following partial rates (as usual, the decay widths of Eqs. (57) are considered 
in units of the free A decay width): 

(62) i?„o = (4.7±2.1)fm3, 

(63) i?pO - (7.9™) fm^ 

(64) i?„i = (10.3 ± 8.6) fm^ 

(65) i2pi = (9.8 ± 5.5) fm^ 

i?„(iHc) = i (i?„o + 3i?„i) = (8.9 ± 6.5) fm^ 
:Rp(iHe) = J {RpO + SRpi) = (9.3 ± 5.8) fm^. 

The errors have been obtained with the standard method, i.e., by treating the data as 

independent and uncorrelated. 

For the ratios of Eq. (60) we have then: 

(66) = 0.6tli, 



(67) ^ = 1.0+1-^ 

Rpl 



1.0- 



while the ratios of the spin-triplet to the spin singlet interaction rates are: 
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The large uncertainties do not allow to draw definite conclusions about the possible 
violation of the A7 =1/2 rule and the spin-dependence of the transition rates. Eqs. (66) 
and (67) arc still compatible with Eq. (59), namely with the AI = 1/2 rule, although the 
central value in Eq. (66) is more in agreement cither with a pure AI = 3/2 transition 
(r ~ 0) or with r ~ 2 [see Eq. (60)]. Actually, Eq. (66) is compatible with r in the range 
— 1/4 40, while the ratio A of Eqs. (60) and (61) is completely undetermined. 
By using the results of Eqs. (62)-(65) one predicts: 

7+11 

^•'-0.7' 
9 O+5.0 

0.95 ±0.92, 

and, by using p2 = 0.001 fm'^ [73], 

rNM(AH) = 0.007 ±0.006. 

The ratio obtained for ^Hc is in agreement with the data of Table XI and with the recent 
determinations discussed in 5'2.2. An accurate measurement of rNM(AH) and F^/Fp for 

and would then provide a test of the weak decay model of Eqs. (57) if the rates 
of Eqs. (62)-(65) could be extracted with less uncertainty from data. 

The compatibility of the data with the AI = 1/2 rule can be discussed in a different 
way: by assuming this rule, we fix Rno/Rpo = 2. Then, by using the observables: 

p 

rNM(AHe), rNM(AHe), :p-(AHe), 

the extracted partial rates are {Rno, Rni, Rn and Rp are unchanged with respect to the 
above derivation): 

i?„o = (4.7±2.1)fm3, 
i?pO = i?«o/2= (2.3±1.0)fm^ 
Rnl = (10.3 ± 8.6) fm^ 
Rpl = (11.7 ± 2.4) fm^ 

These values are compatible with the ones in Eqs. (62)-(65). For pure AI = 1/2 transi- 
tions the spin-triplet interactions seem to dominate over the spin-singlet ones: 



^(iH) = 

^(tH) = 
1 p 

^(iHe) = 
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Moreover, since: 



Rpl 



= 0.9 ±0.8, 



from Eq. (60) one obtains the following estimate for the ratio between the A7 =1/2 
amplitudes: 



{If = {)\\A,,2\\I, = 1/2) 



{If = 1\\A,/,\\I, = 1/2) 



1 

77^-^2.3. 

3.7 



The other independent observables which have not been utilized are then predicted to 
be: 



TNMdH) =0.17 ±0.11, 



(68) ^(iHe) =0.95 ±0.72, 

Ip 

in agreement with the data of Tabic XI, with a for one degree of freedom of 0.31 
(corresponding to a 0.56 ct deviation). This means that these data are consistent with 
the hypothesis of validity of the A7 =1/2 rule at the level of 60%. In other words, the 
A/ =1/2 rule is excluded at the 40% confidence level. On the contrary, if, according 
to the analysis of subsection 5'2.2, one compares with the datum r„/rp = 0.40 ± 0.11 
for /Jle, the A7 =1/2 rule turns out to be completely excluded by the central value of 
Eq (68). 

The observables for which experimental data are not available at present are predicted 
to be: 

^(iH) = 1.3±0.6, 
J- p 

^(1H) = 7.6±6.2, 
J-p 

and, for p2 = 0.001 fm~^, 

rNM(AH) = 0.005 ±0.003. 

We note that the central value of Tn/Tp for in the analysis which enforces the 
A7 = 1/2 rule is considerably larger than the central value obtained in the general 
analysis previously discussed. Thus, the future measurements [68, 69] of this quantity 
will represent an important test of the A7 = 1/2 rule. 
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7. — Non— mesonic decay of polcirized A— hypernuclei: the cisymmetry puzzle 

Lambda hypcrnuclcar states can be produced with a sizeable amount of polarization 
[75] . The development of angular distribution measurements of decay particles (photons, 
pions and protons) from polarized hypernuclei is of crucial importance in order to extract 
new information on hypernuclear production, structure and decay. 

Despite the recent progress discussed in section 5, the reaction mechanism for the 
non-mesonic weak decay does not seem to be fully understood. Indeed, a new intriguing 
problem, of more recent origin, is open: it concerns a strong disagreement between 
theory and experiment on the asymmetry of the angular emission of non-mesonic decay 
protons from polarized hypernuclei. This asymmetry is due to the interference between 
parity-violating and parity-conserving Ap — > np transition amplitudes [76]. The non- 
mesonic rates r„ and Tp are dominated by parity-conserving amplitudes. The study 
of the asymmetric emission of protons from polarized hypernuclei, with its information 
on the spin-parity structure of the Ap np process, is thus supposed to provide new 
constraints on the dynamics of the non-mesonic decay. 

Thanks to the large momentum transfer involved, the n(7r"'", K^)A reaction has been 
used [77, 78], at p^ = 1.05 GeV and small if + laboratory scattering angles (2° < 0k ^ 
15°), to produce hypernuclear states with a substantial amount of spin-polarization pref- 
erentially aligned along the axis normal to the reaction plane. The origin of hypernuclear 
polarization is twofold [75]. It is known that the distortions (absorptions) of the initial 
(7r+) and final (K^) meson-waves produce a small polarization of the hypernuclear or- 
bital angular momentum up to laboratory scattering angles 6k — 15° (at larger scattering 
angles, the orbital polarization increases with a negative sign). At small but non-zero 
angles, the main source of polarization is due to an appreciable spin-flip term in the 
elementary reaction 7r+n AK~^, which interferes with the spin-nonflip amplitude. In 
a typical experimental situation with = 1.05 GeV and 6k — 15°, the polarization of 
the hyperon spin in the free Tr+n — » AK~^ process is about 0.75. 

7T. Spin-polarizaMon observables. - The intensity of protons emitted in Ap np 
decays along a direction forming an angle with the polarization axis is given by (see 
Ref. [79] for more details): 



(69) 



7(e,J) = 7o(J) [1 + ^(6, J)] 



where 



(70) 



Io{J) = 



2J+1 



is the (isotropic) intensity for an unpolarized hypernucleus. In Eq. (70): 



(71 



{J,M) = J2\{F\M\I;J,M)f 



F 
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is the intensity of protons emitted along the quantization axis for a projection M of the 
hypernuclear total spin J. The proton asymmetry parameter, A, can be written in the 
following form [79]: 

(71) A{e, J) = Py{J) Ay (J) cos e. 

The quantity: 



Ay{j) 



3 Y.mMct{J,M) 



which is a property of the hypernuclear non mcsonic decay only, is usually referred to 
as the hypernuclear asymmetry parameter. The hypernuclear polarization Py depends 
both on the kinematics {p„ and Ok) and dynamics of the production reaction. 

In the shell model weak-coupling scheme, Py is directly related to the polarization 
p\ of the A spin in the hypernucleus as follows: 

(72) ..wJ-TTT''.''" '"-■'--i 

[Py{j) if J =Jc + i, 

Jc being the total spin of the nuclear core. It is useful to introduce an intrinsic lambda 
asymmetry parameter a a, which is characteristic of the elementary process Ap np and 
should be independent of the hypernucleus, such that: 

(73) ^(e, J) = pa{J) OA cos e. 
From Eqs. (71) and (72) it follows then: 

aA=< -^^^(-^^ ^J-Jc-k 



Ay{j) if J =Jc + i 



2" 

7'2. Experiment versus theory. - Nucleon FSI acting after the non-mesonic weak 

decay are expected to modify the weak decay intensity of Eq. (69). Experimentally one 
has access to a proton intensity J) which is generally assumed to have the same 

9-dependence as 7(9, J): 

(74) jM(e, J) = I^\J) [1+pa{J) a^{J) cose] . 

The observable asymmetry, ajy (J), which is expected to depend on the hypernucleus, is 
then determined as: 

M.,._ 1 /^(0°, J) -1^(180°, J) 
^ ' ' pa(J)/^(0°, J) + /M(180°, J)" 
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Table XII. - Theoretical and experimental determinations of the asymmetry parameters 
(oa and , respectively). The predictions for have been obtained with different OME 
weak transition potentials and with the direct quark mechanism (DQ). 

Rcf. and Model ^He 

Sasaki et al. [21] 

■K + K + T)Q, -0.68 
Parreno et al. [17] 

TT + p + K + K*+io + rj -0.68 -0.73 

Itonaga et al. [81] 

'K + K + 2'K/p + 2'K/a + uj -0.33 

Barbero et al. [82] 

■K^-p + K + K*+uj + ri -0.54 

KEK-E160 [77] -0.9 ± 0.3 * 

KEK-E278 [78] 0.24 ± 0.22 

KEK E462 (preliminary) [80] 0.07 ± 0.08 

KEK-E508 (preliminary) [80] -0.44 ± 0.32 

* This result correspond to the weighted average (discussed on pag. 95 of Ref. [5]) among 
different p-shell hypernuclear data. 



Until now, four KEK experiments measured the proton asymmetric emission from 
polarized A-hypernuclei. The experiment KEK-E160 [77], which studied p-shell hy- 
pernuclei, suffered from large uncertainties: only poor statistics and energy resolution 
could be used; moreover, the values of the A polarization pA needed to determine the 
asymmetry ojf , had to be evaluated theoretically. More recently, was measured by 
KEK-E278 [78] from the decay of ^He. The values of pA used to obtain were deter- 
mined by observing the asymmetry, = pxa\ , in the emission of negative pions in 
the ^He mesonic decay, after assuming a\ to be equal to the value for the free A — > 'K~p 
decay, a^- = —0.642 ± 0.013. Unfortunately, the small branching ratio and expected 
asymmetry for the mesonic decay of p-shell hypernuclei makes a similar measure- 
ment of p\ very difficult for these systems; even the recent and more accurate experiment 
KEK-E508 [80] had to resort to theoretical estimates for the A polarization in ^^C and 
]y^B. In the other recent experiment KEK-E462 [80], was measured again for aHc, 
but with improved statistics. 

In Table XII we report the results for obtained by the above mentioned exper- 
iments together with available theoretical estimates for oa- While theory predicts a 
negative intrinsic A asymmetry, with a moderate dependence on the hypernucleus, the 
measurements seem to favor positive values for a^ijjie) and negative values for a^Q^^)- 

7'3. Recent developm,ents . - Concerning the above comparison between theory and 
experiment, it is important to stress that, while one predicts aA(AHe) ~ aA(A^C), there 
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Table XIII. - Results of Ref. [64] for the proton intensities [Eqs. (69) and (74)] from the 
non-mesonic weak decay o/^He and 





iHe 

T-M 


M 

"a 


A ^ 
tM 
^0 


M 

"a 


Without FSI 


Jo = 0.69 


OA = -0.68 


lo = 0.75 


OA = -0.73 


FSI and T^^ = 


1.27 


-0.30 


2.78 


-0.16 


FSI and T];^ = 30 MeV 


0.77 


-0.46 


1.05 


-0.37 


FSI and T];^ = 50 MeV 


0.59 


-0.52 


0.65 


-0.51 


FSI and Tp*'' = 70 McV 


0.39 


-0.55 


0.38 


-0.65 


KEK-E462 (preliminary) [80] 
KEK-E508 (preliminary) [80] 




0.07 ±0.08 




-0.44 ±0.32 



is no known reason to expect this approximate equality to be valid for a^. Indeed, 
the relationship between 7(6, J) of Eq. (69) and I^{Q,J) of Eq. (74) can be strongly 
affected by FSI of the emitted protons: this fact prevents establishing a direct relation 
between ua and and to make a direct comparison among results for these quantities. 
In order to overcome this obstacle, an evaluation of the effects of the nucleon FSI on the 
non-mesonic weak decay of iHe, i^B and has been performed very recently [64]. We 
summarize here some results of this investigation, which is the first one evaluating a^. 

The simulated proton intensities turned out to be well fitted by Eq. (74), then one 
can actually evaluate through Eq. (75). In Table XIII we show predictions of the 
OMEf model (a OME model built with the NSC97f potential) for the weak decay and 
observable proton intensities, 1(6, J) and I^{Q,J), respectively. As a result of the 
nucleon rescattering in the nucleus, |aA| > \a^\ for any value of the proton kinetic 
energy threshold: when T^^ = 0, aA/af ~ 2 for iHe and aA/af ~ 4 for fC; \af\ 
increases with T^^ and oa/o^ — 1 for = 70 MeV in both cases. Asymmetries 
rather independent of the hypernucleus are obtained for T^^ ^ 50 and 70 MeV. 
The KEK data quoted in the table refer to a varying between 30 and 50 MeV: the 
corresponding predictions of Ref. [64] agree with the i^C datum but are inconsistent with 
the observation for iHe. 

In conclusion, nucleon FSI turn out to be an important ingredient also when studying 
the non-mesonic weak decay of polarized hypernuclei, but they cannot explain the present 
asymmetry data. Further investigations are then required to clarify the issue. On the 
theoretical side there seems to be no reaction mechanism which may be responsible for 
positive or vanishing asymmetry values. On the experimental side the present anomalous 
discrepancy between different data needs to be resolved. Future experimental studies 
of the inverse reaction pn pA should also be encouraged, since they could help in 
disentangling the puzzling situation. Indeed, the weak production of the A-hyperon 
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could give a richer and cleaner (with respect to the non-mesonic hypernuclear decay) 
piece of information on the lambda-nucleon weak interaction and especially on the A 
polarization observables [83]. 

8. — Summary and perspectives 

In these Lectures we have discussed several aspects of the weak decay of A hypcrnuclci. 
Beyond the mesonic channel, which is observed also for a free A, the hypernuclear decay 
proceeds through non-mesonic processes, mainly induced by the interaction of the A 
with one nucleon or with a pair of correlated nucleons. This channel is the dominant one 
in medium-heavy hypernuclei, where the Pauli principle strongly suppresses the mesonic 
decay. 

Various models have been proposed to describe the mesonic and non-mesonic decay 
rates as well as the asymmetry parameters in the decay of the A hyperon in nuclei. The 
results obtained within these models have been thoroughly discussed. The mesonic rates 
have been reproduced quite well by calculations performed in different frameworks. Also 
the non-mesonic rates have been considered within a variety of phcnomcnological and 
microscopic models. In this context, particular interest has been devoted to the ratio 
Tn/Tp. Indeed, in spite of the fact that several calculations were able to reproduce, 
already at the OPE level, the total non-mesonic width, Tj^m = r„ + rp(+r2), the val- 
ues therewith obtained for Fn/Fp revealed a strong disagreement with the experimental 
data. Although some of these calculations represented an improvement of the situation, 
further efforts were required in order to approach a solution the F^/Fp puzzle. Prom the 
experimental side, recent experiments measured nucleon-nucleon coincidence observables 
with good statistics. Recent analyses of these data, complemented with the theoretical 
estimate of final state interactions, allowed the determination of Fn/Pp values in agree- 
ment with the theoretical expectations. Yet, good statistics coincidence measurements 
of nn and np emitted pairs are further required. These correlation measurements will 
also allow to establish the first constraints on the two-nucleon induced decay width, thus 
testing the various models which have been proposed for the evaluation of the different 
decay channels. 

The analysis of the A J = 1/2 rule in the non-mesonic decay of light hypernuclei 
appears to be feasible, also within a relatively simple phenomenological model: again, 
the measurements of a few delicate, missing decay rates, would help in understanding 
the role of this empirical selection rule. 

As far as the decay of polarized hypernuclei is concerned, the situation is even more 
puzzling. While theory predicts negative values for both the intrinsic asymmetry oa and 
the observable asymmetry a^, with a moderate dependence on the hypernucleus, experi- 
ments seem to favour negative values for aj^(j^G) but small, positive values for ajy (^He). 
Further investigations are then required to clarify the issue. Theoretically, there seems to 
be no reaction mechanism which may be responsible for positive or vanishing ajy values. 
Improved experiments, establishing with certainty the sign and magnitude of for s- 
and p-shell hypernuclei, are strongly awaited. 
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We conclude this work by reminding the reader that hypernuclear physics is 52 years 
old, yet a lot of efforts remain to be done, both experimentally and theoretically, in order 
to fully understand the hyperon dynamics and decay inside the nuclear medium. The 
impressive progress experienced in the last few years is promising and wc hope that it 
deserves a definite answer to the intriguing open questions which we have illustrated 
here. 

* * * 

We warmly acknowledge our colleagues A. De Pace, R. Cenni, A. Parreno and A. 
Ramos, who collaborated with us in obtaining some of the results discussed in these 
Lectures. 
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